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ABSTRACT 

We present soft (0.5-2 keV) X-ray luminosity functions (XLFs) in the Great Observatories Origins 
Deep Survey (GOODS) fields, derived for galaxies at z ~ 0.25 and 0.75. SED fitting was used 
to estimate photometric redshifts and separate galaxy types, resulting in a sample of 40 early-type 
galaxies and 46 late-type galaxies. We estimate k-corrections for both the X-ray/optical and X- 
ray/NIR flux ratios, which facilitates the separation of AGN from the normal/star burst galaxies. We 
fit the XLFs with a power-law model using both traditional and Markov-Chain Monte Carlo (MCMC) 
procedures. A key advantage of the MCMC approach is that it explicitly takes into account upper- 
limits and allows errors on "derived" quantities such as luminosity densities to be computed directly 
(i.e., without potentially questionable assumptions concerning the propagation of errors). The slopes 
of the early-type galaxy XLFs tend to be slightly flatter than the late-type galaxy XLFs although the 
effect is significant at only the 90% and 97% levels for z~ 0.25 and 0.75. The XLFs differ between 
z < 0.5 and z > 0.5, at > 99% significance levels for early-type, late-type and all (early and late- 
type) galaxies. We also fit Schcchter and log-normal models to the XLFs, fitting the low and high 
redshift XLFs for a given sample simultaneously assuming only pure luminosity evolution. In the 
case of log-normal fits, the results of MCMC fitting of the local FIR luminosity function were used 
as priors for the faint and bright-end slopes (similar to "fixing" these parameters at the FIR values 
except here the FIR uncertainty is included). The best-fit values of the change in logL* with redshift 
were A log L* = 0.23 ± 0.16 dex (for early-type galaxies) and 0.34 ± 0.12 dex (for late- type galaxies), 
corresponding to (1 + z) 1,6 and (1 + z) 2 - 3 . These results were insensitive to whether the Schechter or 
log-normal function was adopted. 

Subject headings: galaxies, cosmology, star formation, surveys, x-rays 
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Recently (|Norman et al.l |2004 hereafter N04) pre- 
sented the first X-ray luminosity functions (XLFs) of 
normal/starburst galaxies at cosmologically- interesting 
redshifts (at z ~ 0.3 and 0.7). These XLFs were de- 
rived from the Chandra Deep Field North (2 Ms) and 
South (1 Ms) surveys (hereafter the "CDF" fields). It 
was found that the normal galaxy XLF was consistent in 
normalization and shape to the far-infrared LF, assum- 
ing pure luminosity evolution of ~ (1 + z) 3 . While the 
errors were large due to the limited numbers of galaxies 
in each luminosity bin, the estimated star-formation rate 
(SFR) derived from the XLFs was found to be consistent 
with other SFR measures such as the Ha luminosity. As 
discussed in N04, an interesting aspect of the XLF is that 
the X-ray emission from star-forming galaxies often has 
a large component due to X-ray binaries, particularly at 
hard X-ray energies above 2 keV. Therefore, the XLF is 
in part a probe of the binary stellar mass functions, and 
hence indirectly a measure of the SFR. It is quite possible 
that binaries play a critical role in the evolution of the 
majority of stellar systems, and X-ray emission provides 
one of very few direct probes of such phenomena. 
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The G reat Observatories O rigins Deep Survey 
(GOODS; iGiavalisco et all [2003) is a multi-wavelength 
survey of a sub-area of the CDF fields. This survey 
entails deep imaging by three of NASA's Great Obser- 
vatories: the Hubble Space Telescope, the Chandra X- 
ray Observatory, and the Spitzer Space Telescope, as 
well as extensive photometric and spectroscopic obser- 
vation by ground based facilities. These data allow the 
extension of the N04 results in several new directions. 
First, improved redshift determinations are now avail- 
able for many of the sources. Second, the multi-band 
data have been used to model the spectral-energy distri- 
butions (SEDs) of the sources and estimate their spectral 
types. Therefore XLFs can now be generated as a func- 
tion of galaxy spectral type, which is the primary goal 
of this paper. Finally, we have improved the selection of 
AGN vs. normal/starburst galaxies over N04 by applying 
k-corrections to both X-ray/optical and X-ray/NIR flux 
ratios. Such k-corrections are criti cal to correctly se pa- 
rating AGN from normal galaxies (|Bauer !TaT]l200l . 

A motivation for deriving galaxy type-selected XLFs 
is the fact that there are multiple contributors to the 
X-ray emission of galaxies, namely low-mass X-ray bi- 
naries (LMXRB), high-mass X-ray binaries (HMXRB), 
hot ISM, AGN, and to a lesser extent, individual su- 
pernovae (SN) and massive stars, with the relative con- 
tribution of these exp ected to be dependent on galaxy 
type (sec Ptak 2001, for a review). Specifically, the X- 
ray emission of early-type galaxies is known to be dom- 
inated by LMXRB emission, and in the case of mas- 
sive, gas-rich ellipticals, hot ISM (at the virial temper- 
ature of the galaxy). Late-type spiral galaxies should 
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have significant contributions from both LMXRB and 
HMXRB populations, with the former being associated 
with the older (t > 10 8-9 year) populations. Late-type 
spiral galaxies often exhibit hot ISM due to heating as- 
sociated with recent episodes of star formation. Finally, 
starburst galaxies should have the largest contribution 
from hot ISM (including potentially a superwind out- 
flow) and HMXRB (see iPersic fc Rephaelil 12001 . We 
would therefore naively expect somewhat different evo- 
lution in the X-ray luminosity density of these various 
galaxy types, with the LMXRB contribution following 
the global SFR history of the Universe with a delay of 
the order of the e volutionary time scale of low-mass stars, 
i.e., ~ 10 9 years (|Ghosh fc Whitdl200lh and the HMXRB 
and hot ISM contribution tracking the SFR history in- 
stantaneously (relative to a Hubble time). 

In this paper we are primarily concerned with sources 
with Lx < 10 42 erg s _1 , and hence any AGN present 
would be a low-luminosity AGN (LLAGN). LLAGN are 
found in all galaxy types, with LINER s preferentially 
being found in early- type galaxies (see iHo et al.ll2003l . 
and references therein). As in N04, we assess the im- 
pact of AGN by classifying sources based on a Bayesian 
statistical analysis. In the Appendix we also separately 
address the properties of the X-ray flux ratios (i.e., ra- 
tio of X-ray to optical and near-IR flux) which provide 
important criteria for source classification. Ascertaining 
the presence of any luminosity or rcdshift dependence 
in the X-ray/optical or X-ray/NIR flux ratios would 
also indicate possible evolution of the relative contribu- 
tions of different sources of X-ray flux in galaxies, i.e., 
LMXRB, HMXRB, hot gas, and AGN (for examples of 
stu dy of the evolution of X-ray /optical flux ratio s, please 
see iPtak et aLlfeOOH iHornschemeier et al.ll200l . 

This paper is organized as follows. In Section ?? we 
describe our sample selection and data analysis. The 
results of our analysis are given in Section ??, and we 
discuss the results in Section ??. In the Appendix we 
give the details of our galaxy/ AGN source classification 
procedure, basic statistics concerning the sample, and 
a discussion of the X-ray/optical and X-ray/NIR flux 
ratios. We assume the WMAP cosmology of Hq = 70 
km s^ 1 Mpc" 1 , n rn = 0.3 and fl A = 0.7. 

2. METHODOLOGY 

2.1. X-ray Sample Selection and Redshifts 

The X-ray da t a used in our analysis are taken from 
lAlexander et al.l ((2003) , where the positions and X-ray 
fluxes for sources in both the CDF-N and CDF-S are 
tabulated. Bauer et al. (2007) have carried out de- 
tailed matching of th e GOODS ACS data with the 
lAlexander et al.l ((2003) X-ray catalog, assigning match- 
ing probabilities based on optical faintness of potential 
counterparts, resulting in 263 (185) sources in the CDF- 
N (CDF-S) having a single ACS counterpart. We found 
44 (28) CDF-N (CDF-S) sources with no optical coun- 
terpart, while 3 (4) with multiple counterparts were not 
included in our analysis. We also included the off-nuclear 
X-ra y sources (14 in the CDF- N and 6 in the CDF-S; 
e.g., IHornschemeier et al"1l2004j ). We then used the ACS 
coordinates in the Bauer catalog to match with the pub- 
licly available GOODS spectroscopic redshift catalogs 6 

6 http:/ /www. stsci.edu/science/goods 



and iMobasher et al.l (|2004h photometric redshift cata- 
log to obtain both redshifts and spectral types for all 
sources. Stars were excluded. The main parameters for 
each galaxy are the s oft-band (0.5-2.0 keV) flux (from 
lAlexander et alJ 120031 ). the hardness ratio (defined as 
where H is the 2-8 keV vignetting-corrected count 
rate and S is the 0.5-2.0 keV vignetting-corrected count 
rate; from Alexander et al. 2003), the redshifts, and 
the optical spectral types. We broadly divided the spec- 
tral types into the groups early-type, late-type, and star- 
burst/irregular galaxies, however we only extract XLFs 
for early-type and late-type galaxies since the numbers 
of irregular galaxies were very small. 

2.1.1. Photometric Redshifts 

The photometric redshifts for GOODS fields are esti- 
mated using template fitting technique (Mobasher et al 
1996; Chen et al 1999; Arnouts et al., 1999; Benitez 2001; 
Bolzonella et al 2000, Dahlen et al. 2007). The rest- 
frame Spectral Energy Distribution (SED) for galaxies 
of different types are convolved with filter response func- 
tions of the filters used in photometric observations of 
galaxies. The convolved SEDs, shifted in redshift space, 
were then fitted to observed SEDs of individual galaxies 
by minimizing the % 2 function 

X 2 = ^UW obs -uFt emplate )lo-*)] 2 

where the summation, i, is over the passbands (i.e. num- 
ber of photometric points) and n is the total number of 
passbands. F£ bs and F^ emplate are, respectively, the ob- 
served and template fluxes at any given passband. a % is 
the uncertainty in the observed flux and a is the normal- 
ization. The redshift and SED (i.e. spectral type) cor- 
responding to the minimum \ 2 value for a given galaxy 
were then assigned to that galaxy. We used priors based 
on luminosity functions (LFs). The main effect of a LF 
prior is to discriminate between cases in which the red- 
shift probability distribution function, which identifies 
the most likely redshift, shows two or multiple peaks 
(i.e. more than a single optimum redshift) due to con- 
fusion between the Lyman break and the 4000 A break 
features. The absolute magnitudes of the object at the 
redshift peaks can then discriminate between these pos- 
sibilities. Absorption due to intergalactic HI is included 
using the parametrization in Madau (1995). 

We use template spectral energy distributions (SEDs) 
for normal galaxies consisting of E, Sbc, Scd and Im from 
Colman et al (1980) and two starburst templates from 
Kinney et al (1996)-(SB2 and SB3). To increase the spec- 
tral resolution, we construct intermediate- type templates 
by using the weighted mean of the adjacent templates 
and interpolate between them. This is done by defining 
five intermediate-type templates between the main spec- 
tral types used. Therefore, we use a total of 31 SED 
templates in this study. Photometric redshifts were then 
measured using the the observed photometry in GOODS- 
N (UBVRisJK) and GOODS-S (UBVRilHK). 

2.1.2. Spectroscopic Redshifts 

For a subset of the galaxies, we use the available 
spectroscopic redshifts. The majority of the spectro- 
scopic redshifts were derived from Keck DEIMOS data 
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for the CDF-N (see IWirth et al1l2004 and references 
th erein) and VLT FORS2 or VIMOS data fo r the CDF- 
S dSzokolv et al.ll200l IVanzella et al.ll200l . The limit- 
ing magnitude for spectroscopic redshift determination 
was typically Rab ~ 24. If a spectroscopic redshift was 
not available then a photometric redshift was used, with 
the photometric redshift limiting magn itude typically be- 
ing Rab ~ 25 (|Mobasher et al.ir2004h . In cases where 
a quality assessment was available for the spectroscopic 
redshift and was considered to be poor, and a photo- 
metric redshift with error S z < 0.2 was available 7 , the 
photometric redshift was used. In cases where there was 
no spectroscopic redshift quality given and a photomet- 
ric redshift with error 5z < 0.1 was available, the pho- 
tometric redshift was used. This resulted in 204 (129) 
and 157 (116) total (spectroscopic) redshift determina- 
tions for the CDF-N and CDF-S sources, respectively 
(i.e., 73% and 82% of the GOODS X-ray sources with 
a unique optical counterpart have a redshift estimate). 
The redshift distributions are shown in Figure [T] 

In cases where both a photometric and spectroscopic 
redshift were available, the mean absolute redshift differ- 
ence was 0.09 (CDF-N) and 0.13 (CDF-S) after removing 
a small number of outliers with Az = |2; spcc — z p hot > 0.5. 
In Figure [2] we plot the difference between the photo- 
metric and spectroscopic redshift for sources with both 
measurements. The errors plotted in the figure are 
based solely on the photometric redshift errors. This 
figure shows that our typical redshift uncertainty is 
Az ~ 0.1 — 0.2 for z < 1.5, and increases significantly 
for z > 1.5. We conservatively limit our analysis to 
z < 1.2. This leaves a total of 148 sources in the CDF- 
N and 95 sources in the CDF-S regions, of which 138 
and 81 have soft-band X-ray detections. Note that 131 
(CDF-N) and 75 (CDF-S) of the X-ray sources have 
Lx < 1 x 10 43 ergs s" 1 , the highest X-ray luminosity 
considered in the X-ray luminosity functions discussed 
here. 

2.2. K-Corrections 

The soft-band X-ray fluxes, listed in Alexander et al. 
(2003), are based on the observed count rates and a count 
rate to flux conversion computed from an effective photon 
index (based on band ratios). This photon index (T) was 
used to "k-correct" the observed fluxes, i.e., Fx.kcor = 
(1 + z) r - 2 F x ,obs, and L x = ^d\F x .kcor where d L is 
the luminosity distance. In Appendix A we present the 
sources, the adopted redshifts and X-ray luminosities. 

For a given band Q, we define Fq to be vF v at the 
central frequency of the band vq. Since the GOODS 
magnitudes are in the AB system, log(Fx/FQ) is given 
by log Fx + 0.4(Q + 48.6) - \ogVQ, e.g., in the case of 
R band, \og{F x /F R ) = \ogF x + 0AR AB + 4.8. The 
optical magnitudes were k-corrected by interpolating be- 
tween the magnitudes whose central wavelengths bracket 
(1 + z)X, where A is the central wavelength of the fil- 
ter of interest (e.g., 4400A for B-band). The band- 
pass corre ction term (1+z) wa s also included in the k- 
correction ([Blanton et al.ll2003f) . In Figure [3] we show the 
k-corrections applied to R and K s band flux ratios. Note 
that in the case of K s band magnitudes, the interpolation 

7 Sz is defined as the 68% uncertainty on the photometric redshift 
derived from the posterior probability 



discussed above necessarily amounts to an extrapolation 
since K s is the longest wavelength band in the dataset 
and every source considered here has z > 0. The correc- 
tions are plotted separately for each spectral type, where 
it can be seen that the largest k-corrections occur for B 
band fluxes from early- type galaxies, which is not sur- 
prising considering their red color. Also, as mentioned 
in Bauer et al. (2004), the k-corrections to the flux ratios 
are dominated by the optical correction rather than the 
X-ray correction. 

2.3. Luminosity Function Estimation 

We construct bin ned luminosity functions using the es- 
timator described in lPage &: Car rcra (2000). Briefly, this 
assumes negligible variation in both the luminosity func- 
tion and its evolution across each ALAV(z) bin, where 
V(z) is the co-moving volume: 

J C(L,z)V(z)dzdL (1) 

N is the number of sources in the XLF bin bounded by 
Lmim L max , z m in and z max , and C(L,z) is a complete- 
ness correction. Note that z max is a function of luminos- 
ity since the limit is chosen at the highest observable red- 
shift given the limiting flux (or, equivalently, C(L, z) = 
where F(L,z) < Fu m ). The error for each XLF bin is 
accordingly derived from the Poisson error on the num- 
ber of galaxies in each bin (note this assumes negligible 
statistical and systematic error in the volume integral 
terms), and for bins w ith no galaxie s we used N= 1.841 
as the Id upper limit (jGehr cls 1986). Small luminosity 
function bins are preferred, however larger bins minimize 
Poisson noise and reduce the effect of luminosity uncer- 
tainties within a bin (due to redshift uncertainties). For 
example, for the uncertainty in luminosity resulting from 
redshift error to be on the order of the bin size, bin sizes 
of AlogL = 0.5 are required for sources with z > 0.4 
and redshift errors of 0.1 or less. Given the small num- 
ber of sources, we divide our sample into just two redshift 
bins, z < 0.5 and 0.5 < z < 1.2 . We used the soft-b and 
GOODS sky coverage shown in lTreister et alj (|2004h . 

2.4. Completeness Correction 

The completeness correction C(L,z) is given by the 
product of the X-ray detection completeness, the prob- 
ability of a given X-ray source having an optical/NIR 
counterpart, and the probability of a counterpart having 
a redshift. The latter two terms reduce to the probabil- 
ity of an X-ray source having a redshift. For the X-ray 
detection completeness we used the results of simulations 
originally performed for the full CDF areas discussed in 
Bauer et al. (2004), where the effective solid angle for 
each source was computed (i.e., the maximum solid an- 
gle over which the source could have been detected) . The 
ratio of the effective solid angle to the (total) CDF solid 
angle at the flux of the source is then an estimate of 
the completeness for sources at a similar flux and po- 
sition. We then computed the completeness as a func- 
tion of flux by multiplying the mean per-source com- 
pleteness (i.e., effectively averaging over off-axis angle) 
and the fraction of GOODS x-ray sources with a red- 
shift in the given flux range (the bin sizes were 0.5 dex 
for the GOODS-N sources and 0.75 dex for the GOODS- 
S sources), shown in Figure [4j In both the GOODS-N 
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Fig. 1. — The redshift distributions for the CDF-N (left) and CDF-S (right) fields. The spectroscopic redshift distributions are shown in 
blue with dot-dashed lines. In cases where the spectroscopic redshift quality was unknown or flagged as low and the photometric redshift 
error was small (see text for details), the photometric redshift was adopted. This distribution of sources for which a spectroscopic redshift 
was adopted is shown in red with dashed lines. 
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plotted as a function of z spec for CDF-N (left) and CDF-S (right). In the CDF-S plot, the sources with poor 
spectroscopic redshift determinations (often based on a single line identification) are shown in red and marked with crosses) The errors 
plotted for the points are based solely on the photometric redshift error. 



and GOODS-S the completenesses ranged from ~ 60% 
at Eq. 5_2.o kcV = 2 x 10 -17 erg cm -2 s _1 in the North 
and Fo.5-2.0 kcV = 7 x lCT 17 in the South to 100% for 
Fq. 5-2.0 kcV > 10 _14 erg cm~ 2 s -1 , although with errors 
on the order of 20% in each flux bin. We fit the complete- 
ness with a quadratic function, both to avoid the impact 
of statistical fluctuations and to have an analytic expres- 
sion for use in calculating Eq. 1. The best-fitting relation 
was C(F X ) = 14.56 + 1.66 log + 0.0495(logF x ) 2 and 
is shown in Figured] 

2.5. Galaxy Classification 

N04 classified sources as AGN1 (broad-line AGN), 
AGN2 (narrow- line AGN), and "galaxy" (i.e., nor- 
mal/starburst galaxy) using a prior based on the ob- 



served distributions of X-ray luminosity (\ogLx), the X- 
ray-to-optical flux ratio (log Fx+R) , and X-ray hardness 
ratio derived from a subset of the sources where the op- 
tical classification was secure. These distributions were 
used as priors along with the likelihoods of these mea- 
sured values in computing the posterior probabilities. A 
source was then classified as a galaxy if the galaxy prob- 
ability exceeded the AGN1 and AGN2 probabilities. 

Here we employ a similar algorithm, discussed in Ap- 
pendix A . However, there are a number of differences 
with the approach in N04, including use of the X-ray/K- 
band flux ratio, inclusion of k-corrections in the compu- 
tation of X-ray/optical and X-ray /NIR flux ratios, and a 
more conservative galaxy selection threshold. Note that 
our priors are based on the spectroscopic data in Szokoly 
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Fig. 3. — The k-corrections applied to the X-ray flux ratios for the R (left) and K s (right) bands, for both of the GOODS fields. The 
points show Alog[F x / F opt ] = log [(Fx / 'F op t) z =o] — log[(Fx/F op t)^], where log[(Fx/F op t)z=o] is the flux ratio after applying k-corrections 
and log[(Fx/F opt ) 2 ] is the "observed" flux ratio. Early-type sources are marked with (black) diamonds, late-type sources are marked with 
(red) pluses, and irregular /star burst sources are marked with (blue) asterisks. The solid and dotted lines show the X-ray k-correction for 
f 1 = 1.8 and T = 1.0, respectively. 

resulting in HR<0 in AGN2/starburst composites (the 
case for 4 out of the 7 AGN2 with HR<0). Our final 
normal/starburst galaxy sample consisted of 64 CDF-N 
sources and 23 CDF-S sources. We also consider an "op- 
timistic" normal/starburst galaxy sample which includes 
all sources that were not classified as an AGN. This sam- 
ple therefore includes sources with ambiguous classifica- 
tion, similar to the original selection in N04. The opti- 
mistic normal/starburst sample contains 98 CDF-N and 
54 CDF-S sources. 

2.6. XLF Model Fitting 

In order to quantify differences among the various 
XLFs discussed here, we fit the XLFs with several mod- 
els (note that all luminosity functions discussed here are 
binned). The simplest model is a linear function in 
log^-logLx space, i.e., \og<fi(Lx) — alogLx + b. We 
performed linear fits with two conventional approaches, 
namely using the survival analys is linear reg r ession (here- 
after Method 1) discussed in llsobe et all (|1986l) . and 
least-squares fitting (using the QDP program) after ex- 
cluding the upper limits (Method 2). The linear regres- 
sion approach has the disadvantage of not explicitly in- 
cluding any error information (i.e., the algorithm only 
takes as input the value of each detection or limit), and 
the confidence level used in computing an upper limit is 
arbitrary. 

2.6.1. Markov-Chain Monte Carlo Fitting 

Bayesian parameter estimation has been gaining pop- 
ularity in various fiel ds of astronom y, most notably in 
cosmology (see, e.g.. ISpergell l2007h . A key advantage 
of Bayesian parameter estimation is that common sta- 
tistical issues, such as the treatment of upper limits in 
fitting and propagation of errors, are inherently handled 
properly. In addition the probability distributions for 
parameters of interest are computed and can be shown 
rather than simply a summary such as the 68% confi- 
dence interval. This latter point is particularly relevant 
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Fig. 4. — Completeness corrections as a function of flux cal- 
culated for the GOODS-N (crosses) and GOODS-S (diamonds) re- 
gions. The line shows a quadratic fit to the combined completeness 
corrections (i.e., to be applied to either data set). 



et al. (2004) which covers the entire CDF-S area. In 
Szokoly et al. only R and K band values are listed, and 
accordingly k-corrcctions could not be interpolated as 
was done with the GOODS data (see §2.2). However, we 
estimated the spectral type of the sources based on their 
R-K color, and applied a mean k-correction derived for 
each type. 

In the case of X-ray hardness ratio, the means and 
standard deviations of the prior samples were —0.19 ± 
0.46 for galaxies, 0.16 ± 0.37 for AGN2 and -0.51 ± 0.05 
for A GN1, very similar to t hose used in N04. For compar- 
ison, iPeterson et al.l (|2006f ) artificially redshifted a sam- 
ple of local galaxies to z = 0.3 and found that hard- 
ness ratios of AGN2 was split evenly around HR=0. 
They noted that soft emission from starbursts could be 
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when there are relative minima in the fitting statistic, 
in which case the meaning of traditional error bar is not 
well defined. 

The basis of Bayesian parameter estimation is the com- 
putation of the posterior probability distribution: 



\D) 



p(0)p(D\0) 
p(D) 



(2) 



where 9 is the vector of model parameters (e.g., 9 — {a, b} 
for the linear model), D represents the data, p{9) is 
the prior probability distribution for the parameters, 
p(D\9) is the likelihood function, and p{D) is a nor- 
malization constant (i.e., since it does not depend on 
the parameter values). p{D) is given by J p(6)p(D\9)d9. 
In practice, computing p(9\D) is computationally diffi- 
cult since it requires an n-dimensional integral when fit- 
ting an n parameter model. The Markov-chain Monte 
Carlo (MCMC) procedure circumvents this by perform- 
ing a directed random walk thro ugh the parameter space 
van Dvk et al.ll2(i0llForJ20f)6h . Here we assume Gaus- 
sian priors for each parameter (with mean and standard 
deviation pg i and ag i for parameter 9i), and chose very 
large ag i values for the case of a flat (uninformative) 
prior. The likelihood function for the number of counts 
in an XLF bin is the Poisson distribution, giving 

n 

P (9\D) ^Y[G(9 Z \ 

i=l 

m 

xY[Pois[N J \ ( f > {9 1 ,..,9 n ,L Xj ,z 3 )V{L Xj ,z 3 )} (3) 
i=i 

where there are n model parameters and m XLF bins. 
G(9i\nen 00;) is the Gaussian prior for parameter 9i. 
Nj is the number of sources and V(Lxj, Zj) is the co- 
moving volume corresponding to the jth bin (see Eq. 1). 
4>{9\, ..,8 n , Lx.j, zj) is the model XLF evaluated at Lxj 
and Zj. The Pois[) terms give the likelihood of detect- 
ing Nj galaxies in XLF bin j given an expectation of 
<fiV. In addition to the simple linear model we also fit- 
ted Schechter and log-normal functions to the individual 
XLFs. We used the Metropolis-Hastings MCMC algo- 
rithm, in which a "proposal" distribution 8 is used to 
guide the variatio n of the parameters (see iFordl 120051 : 
iRussell eT al. 2007). In this procedure random offsets for 
each parameter are drawn from the proposal distribu- 
tion, and accordingly the step sizes (i.e., the Gaussian 
as) are preferably on the order of the final error for that 
parameter. A step is accepted if the probability of the 
model given the new parameter values is higher, and also 
at random intervals when the probability is lower (i.e., 
occassionally the fit is allowed to proceed "downhill" to 
avoid relative minima). For a given run, three chains 
were produced with a length of at least 2 x 10 6 itera- 
tions, and the parameter step sizes were adjusted during 
the first chain to achieve acceptance rates in the range 
0.35-0.4. Only the last chain in a given run was used 
for analysis. At least 10 runs were performed for a given 
fit and we computed the convergence R statistic from 
iGelman et"aH (|2004D . R values < 1.2 indicate conver- 
gence, and in every case the R values were < 1.01. In 

8 We used the Gaussian distribution as the proposal distribution 
as is common practice. 



practice, the linear model parameters a and b were highly 
correlated, which leads to i neffici ent MCMC fitting. As 
discussed in IGelman et al.l (|2004D . we addressed this by 
instead fitting for ax, a,2 and &, where a — a% + 0,2b. The 
initial values of di and 02 were determined from a linear 
regression fit to the output of a short MCMC run, and 
the value of a<i was held approximately fixed by using a 
tight prior (a\ and 02 are "nuisance" parameters and are 
not discussed further) . The initial parameter values for a 
given run was chosen randomly and the range of the al- 
lowed values was ±2x the initial step size from the prior 
mean. The dispersion between the best-fitting parameter 
values from the runs in a given fit was always < 10% of 
the final error in the parameters, showing that the chains 
converged independently of the starting values. For each 
parameter of interest the chain values were binned into 
normalized histograms, which represent the marginalized 
posterior probability distribution of the parameter. The 
peak value (i.e., the mode) was taken as the best-fit value 
and the 68% confidence interval was derived from the 
probability density by steeping from the peak in the di- 
rection of the smallest decre ase until the integ rated area 
equaled 68%, as discussed in lKraft et alJ ([1991). In order 
to assess how constraining the prior is, we computed the 
ratio of the prior a values to the standard deviation of the 
MCMC parameter values (i.e., the 68% error in the case 
of a Gaussian posterior probability distribution). In ta- 
bles showing the fit results we marked parameter values 
and errors as having a tightly-constraining prior when 
this ratio is <1.1 and as a moderately-constraining prior 
when it is between 1.1 and 2.0. 

2.6.2. Fit Quality 

Unfortunately the MCMC procedure does not directly 
provide a model probability estimate since this requires 
the normalization of Eq. 2 to be computed. Since ef- 
ficient model selection is complicated (and the subject 
of active research, see, e.g., Trotta 2007), we defer this 
and here we compute x 2 for the models after excluding 
upper limits. This also has the advantage that it can be 
consistently applied to all of the model fitting methods 
discussed here in the case of the linear model. 

3. RESULTS 

3.1. Linear Model Fits 

In Figures [5][9] we show the XLFs derived from the 
GOODS fields based on the early-type, late-type and to- 
tal early + late-type galaxy samples along with the lin- 
ear model fits. We also show the fits to the "optimistic" 
galaxy sample and the N04 XLFs for comparison. The 
best-fit parameters and errors are given in Table [1] along 
with the associated x 2 values and the probability with 
which the model can be rejected (again with the caveat 
that this probability does not include upper limit data). 

As is evident from the fit results, the MCMC fitting 
was not very sensitive as to whether the Method 1 or 
Method 2 fit values were used as the prior means. Inter- 
estingly, in every case the MCMC fits resulted in curves 
that were intermediate to the Method 1 and Method 
2 fits, and tended towards the Method 1 results when 
upper- limits were constraining (e.g., the early- type and 
late-type galaxy XLFs) and towards Method 2 when 
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Fig. 5. — Full-sample (early + late-type galaxy) XLFs from the GOODS fields, for the redshift intervals z<0.5 (left) and 0.5 < z < 1.2 
(right). The solid (black), dashed (red), dot-dashed (blue) and long-dashed lines (green) lines respectively show the best-fitting models 
from Method 1 (survival analysis linear regression), Method 2 (least-squares excluding upper-limits), MCMC fitting with Method 1 results 
used to initialize the fit, and MCMC fitting with Method 2 results used to initialize the fit. 
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Fig. 6. — "Optimistic" sample (early + late-type galaxy) XLFs from the GOODS fields, for the redshift intervals z<0.5 (left) and 0.5 < 
; < 1.2 (right). Lines are as in Figure|5] 



the no upper limits were used or the upper limits were 
not constraining (e.g., the Normal et al XLFs). The 
marginalized posterior probability distributions for a and 
6 were close to Gaussian, resulting in nearly symmetric 
errors. In general the linear models cannot be rejected 
at high confidence (i.e., > 3er). Hereafter, the fit results 
presented refer to MCMC fitting. 

3.2. Joint Linear Model Fits 



In Figures [TOlfTTI we show the results of fitting corre- 
sponding pairs of XLFs simultaneously (i.e., the low and 
high redshift XLFs for a given sample, and the early and 
late-type galaxy XLFs at low or high redshift), explicitly 
fitting for the offsets in a and b between the XLFs. Fig- 
ure[T2]shows the early and late- type XLFs similarly being 
fitted jointly. Here also the posterior probabilities were 



mostly Gaussian in shape, and the best-fit parameter val- 
ues and errors are given in Table [5] These fits resulted 
in parameters for the given XLFs that were equivalent to 
those obtained by fitting the XLFs seperately with the 
linear model, however here we are deriving the probabil- 
ity distribution for the offsets in slope and intercept. As 
discussed above, this approach avoids the need to propa- 
gate errors in determining the significance of the change 
in the linear model parameters between XLFs. How- 
ever, in these fits the final errors on Aa and Ab are very 
similar to adding the errors obtained from the individ- 
ual fits in quadrature, which is perhaps not surprising 
given that the posterior probabilities for these parame- 
ters were nearly Gaussian. In every case, the best fit 
values of Aa were < and the best fit values of Ab were 
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The N04 XLFs for z<0.5 (left) and 0.5 < z < 1.2 (right), with linear fits shown as in Figure[5] 
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Fig. 
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z < 0.5 early-type (left) and late-type (right) galaxy XLFs from the GOODS fields. Lines are as in Figure [5] 
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> 0. In order to estimate the significance of a change 
in the XLF, we determined the fraction of simulations 
in which Aa < and Ab > 0. Since in some cases this 
probability is very small, we performed 50 runs with a 
very large chain length (5 x 10 6 ) for these fits. These 
probabilities are also listed in Table [H 

3.3. Log-Normal and Schechter Fits 

We fit the XLFs using the log-normal (Saunders et 
al. 1990) and Schechter (Schechter et al. 1976) func- 
tions since these functional forms fit the FIR (and hence 
star-forming galaxy) and optical luminosity functions of 
galaxies well. Specifically the functional forms were: 



exp 



~i log2(1 ^ 



L 

L* 



o(I) = ln(1<no (i- 



"MM — 



(4) 



(5) 



where in both cases the units of (f> are galaxies Mpc -3 
logL -1 . For both of these sets of fits we placed very 
weak priors on log L* and log 4>* . We first fitted the FIR 
LF published in Saunders et al. (1990) in order to ob- 
tain prior information based on the local star-forming 
galaxy luminosity function, and to validate our method- 
ology. We assumed that the errors listed for each LF 
point were Gaussian and excluded upper limits. Our 
results are shown in Table [3l along with the original re- 
sults of Saunders et al. (1990) and the results of fitting 
the log-normal function to a more recent FIR sample in 
Takeuchi et al. (2001). Our results are consistent with 
the fitting based on traditional methods within the er- 
rors. Since the Takeuchi et al. LF fits were based on 
more recent data we set the log-normal fit prior means 
to their values, with log L* scaled to the X-ray band using 
logLx = logL 60 Mm — 3.65 (Ranalli et al. 1990), and by 
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Fig. 9. — 0.5 < z < 1.2 early-type (left) and late-type (right) galaxy XLFs from the GOODS fields. Lines are as in Figure[5] 
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Fig. 10. — Joint linear fits to the low-z and high-z XLFs for the full (early + late-type) galaxy sample (left) and the optimistic galaxy 
sample (right). The low-z XLF points are marked with (black) squares and the high-z XLF points are marked with (blue) triangles. The 
dashed (black) lines and dot-dashed (blue) lines respectively show the fit to the low-z and high-z XLFs. 



(1+z) 3 where z = 0.25 for the z < 0.5 XLFs and z = 0.75 
for the 0.5 < z < 1.2 XLFs. In the case of the early-type 
or late-type galaxy XLFs we further reduced the prior 
mean for (f>* by a factor of 2 since these samples were 
~ 50% of the total galaxy sample. We conservatively as- 
sumed the MCMC errors in Table [3] which were similar 
to, but larger than, the Takeuchi et al. errors. We set 
the prior standard deviations to 5 times these errors for 
a and a and 50 times these errors for logL* and log</>* 
(effectively weak). In practice, however, different values 
of a were preferred by the data, and so we refit the data 
after setting the prior mean (and initial starting values) 
to the best-fitting a values from the original fits. The 
fits were not stable if we allowed the slope prior widths 
to be broader, and in the case of the early-type galaxy 
XLFs, the prior for a was required to be 50% smaller in 



order to result in stable fits. 

The log-normal XLF fits to the XLFs are shown in 
Figures [T4]fT8l with the best-fit parameters and errors 
given in Table |H The posterior probability distributions 
for the fit parameters for the z < 0.5 early-type sam- 
ple XLF are shown in Figure [THl and the probability 
distributions for the fit parameters from the other log- 
normal fits have similar shapes. Note that using a tighter 
prior for a for the early-type galaxy fits did not change 
the best-fit parameter values significantly but resulted in 
~ 25% smaller errors for log L* and log </>* . The posterior 
probability distributions for a were obviously completely 
dominated by the prior distribution. This is, neverthe- 
less, an improvement over simply "fixing" fit parameters, 
which is the common procedure when parameters are not 
sufficiently constrained by the data since the width of the 
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TABLE 1 
XLF Linear Fits 



Sample 


Method 


a 


b 


X 2 /dof 


A 


__ar_y-typc galaxies, low-z 


i 

i 


— u. iu _e u.ii 


9^ O -1- ZL ^ 

ZO.U It _:.0 


Q £o /9 
O. OO / Z 


OO. z /o 




2 


-0 53+ 013 

u - o,:5 -0.13 


18 3+ 54 


6.34/2 


95.6% 




3* 


u - ' -0.09 


25 4 + 5'2 
iJl -3.3 


3.66/2 


83.9% 




_ t 


-0.72I .i 


25.6±|;| 


3.70/2 


84.2% 


T 

Late- type galaxies, low-z 


i 

i 


n 00 _|_ n 1 /i 

— u.yy it U.14 


OO. O it 0.0 


Q 7S /9 
0. to/ Z 


O4.O/0 




2 


-0 7Q+ 011 


28 6+ 4 7 

z °-°-4.6 


5.55/2 


93.6% 




3* 


-0 91+ 09 
"•"-0.11 


33 6 +4 1 

oo.u_ 3 g 


2.95/2 


77.0% 




■■> f 


U - Mi -0.08 


33 - 6 -4.0 


2.97/2 


77.2% 


All galaxies, low-z 


i 

i 


n on -t- n 1 1 


^9 ^ -1- ZL R 
OO.O It -i.D 


o.yo / z 


yo.D /o 




2 


-0 67+ 09 
u '°'-0.08 


24 1+ 3 ' 4 


12.47/2 


99.5% 




3* 


-0 81+ 06 


29 a+3-3 


6.64/2 


96.2% 




3 1 


-0-82-0.08 


30.1±fl, 


6.64/2 


96.2% 


All 

All galaxies (optimistic), low-z 


i 

i 


n ni 4. n 1 1 
— u.yi ± u.11 


9Q Q _l_ yl £? 

00.0 it 4.0 


1 OA /Q 

lo.yo/o 


00 o*% 
yy.y /o 




2 


-0 67+ - 07 
u -°' -0.07 


24 4+ 2 9 


19 03/3 


>99 0% 






n -7+0.06 
u - ' 1 -0.06 


28 1+ 2 ' 4 

io - -2.4 


9.85/3 


98.0% 




31 


„ ~~+0.06 

-0-77+^;^ 


28.2+2's 


9.70/3 


97.9% 


iMorrnan ex ai. ^_,uu_:j, iow-z 


i 

i 


O Q9 -1- O 07 
— u.yz It U.U 1 


9-1-9 8 


90 OA I A 


^QQ 0% 

^yy.y /o 


2 


-0 76+ - 05 
u - '"-0.05 


97 c+2.2 
z 1 -0_2.i 


13.85 /4 


99.2% 






n 77 +0.05 
u - ' ' -0.04 


28 2+ 2 ' 
z °- z -1.8 


10.56/4 


96.8% 




■-> + 

31 


„ „ o +0.04 

-0-78to.o6 


no ,1 + 1.7 

28.4+2.1 


10.43/4 


96.6% 


__ar_y-typc galaxies, lii-z 




1 -1- 1 


K1 K _|_ /I O 


1 17/9 


AA 




2 


1 91+0.21 
J - a — 0.21 




1.25/2 


46.3% 






-1 28+ ' 20 


49 0+ 9 ' 1 


1.07/2 


41.4% 




3 1 


1 n^+0.18 

-J- 27 I .23 


ao -1+9.5 
48 - 4 -7.5 


1.07/2 


41.3% 


Late- type galaxies, hi-z 


1 


1 an -1-0 10 
— i.yu it u.iu 


7 A ^ -1- A 9 

f _:.0 III 


U. OO / 1 


oo.y /o 




2 


-1 79+ 029 

x - ' a -0.29 


70 1+ 11 - 7 
,u - 1 -11.7 


0.52/1 


51.5% 






-1 85 +a29 
l.OO_ 23 


72 6+ 12 - 4 
,z -°-8.8 


0.54/1 


52.3% 




3 1 


_ s .-,-)_0.19 
l.Od_ 32 


rr. -+13.0 

71.7_g { 


0.49/1 


50.5% 


All galaxies, hi-z 


1 


— 1 .OO IE U.UO 


(\A (\ -+ *? 1 
o*±.o zn 0.0 


QQ /9 

u.yy / z 


oo.y /o 




2 


± - J± -0.16 


c-o Q+6.7 
do - a -6.7 


0.62/2 


26.8% 




3* 


-L . OO n 1 1 

— U. 14 


61.7+H 
— 7.3 


0.64/2 


27.2% 




3t 


1 r l9 +0.13 


59.3l|i 


0.60/2 


25.9% 


All galaxies (optimistic), hi-z 


1 


-1.50 ±0.17 


58.7 ±6.9 


11.75/2 


99.4% 




2 


_1 94+012 
0.12 


47.8+gi 


13.83/2 


99.6% 




3* 


, oq+0.11 

i -00_g 09 


51.8t 4 4 '? 


9.16/2 


98.7% 




3t 


_i 34+0.10 
l-^-O.lO 


52.21J 9 


9.12/2 


98.7% 


Norman et al. (2004), hi-z 


1 


-1.42 ±0.10 


55.1 ±4.2 


15.25/3 


99.8% 


2 


1 1 7 +0.09 
± - 1 '-0.09 


45.0+f;; 9 


4.10/3 


74.9% 




3* 


-1 18+ - 10 


45.113- 6 


3.76/3 


71.1% 




3t 


1 10+O.O8 
-J.-J-8_0.10 


45.213-6 


3.76/3 


71.2% 



Note. — Fit parameters are given for log<f> = alogLx + b. Method 1: linear regression, 2: fit to data excluding upper-limits, 3*: 
MCMC fit using Method 1 results to initialize the fit, 3^: MCMC fit using Method 2 results to initialize the fit (see text). p x 2 gives the 
X 2 probability at which the model fit can be rejected (note that x 2 is computed excluding upper limits). 



priors are scaled from the "physical prior" of the 60 /jm 
LF fit results. A broad tail in the probability densities 
toward low values of log (j)* and high values of log L* is 
due to degeneracy between these parameters 9 . This also 
can be seen in Figure [20] where the MCMC draws are 
plotted. 

We computed the luminosity density, p = 
J 4>(L) Ld log L, for each MCMC draw by numeri- 
cally integrating the log-normal function over the range 
10 37 < Lx < 10 43 erg s -1 . The posterior probabilities 
for p are shown in Figure [21] Note that this results in a 
statistically-correct estimation of p and its error since no 
(usually questionable) propagation of errors is required. 

9 Note that the prior is not biasing this result since in both cases 
the prior peaks on the opposite side of the probability density to 
the broad tail. 



In the case of the Schechter function fits we assumed 
a = 1.0 with a prior width of 50% (i.e., prior p a — 

1.0, cr a = 0.5) since the GOODS-S J-band 0.1 < z < 
0.5 luminosity function a values where in the range 
-1.4 < a < -0.5 (Dahlen et al. 2005). We also used 
the early-type galaxy 0.1 < z < 0.5 J-band LF fit pa- 
rameters from Dahlen et al. (2005) for MJ (-22.97) and 
(/)* (8.6 x 10~ 4 Mpc~ 3 mag^ 1 ) to estimate the initial 
XLF fit parameters (and prior mean values) for log L* 
and (f>* . We converted Mj to the X-ray band by com- 
puting the mean k-corrected value of log Fx + 0.4 J + 5.1 
(see also Appendix C) to be -3.3 for normal galaxies. (f>* 
was rescaled by 2.5 to be in the units of Mpc~ 3 dex -1 . 
As with the log-normal fits the prior mean for <p* was 
reduced by a factor of 2 for the early-type and late-type 
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Fig. 12. — Joint linear fits to the low-z early and late-type XLFs (left) high-z early and late-type XLFs (right). The early-type XLF 
points are marked with (black) squares and the late-type XLF points are marked with (blue) triangles. The dashed (black) lines and 
dot-dashed (blue) lines show the fit to the early and late-type XLFs. 

galaxy sample XLFs. The best-fitting Schechter mod- 
els are also shown in Figures [T4lfT8l with the best- fitting 
parameter values and errors given in Table [5] The pos- 
terior probability densities for logL*, log<^*, a and p are 
shown in Figurel22lfor the z < 0.5 early-type galaxy sam- 
ple, and again other Schechter fits resulted in posterior 
probabilities with roughly similar shapes. 

We also fitted the low and high-z XLfs simultaneously, 
in this case only allowing the log L* to vary between the 
low and high-z models. This is, by definition, the case 
of pure luminosity evolution (PLE) . The results of these 
fits are shown in Figures and the fit parameters are 
listed in Tables [5] and [JJ In all cases the posterior proba- 
bility for A log L* was nearly Gaussian, and we show the 
cases of the log-normal fits to the early-type and late- 
type galaxy XLFs in Figure [25] 
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TABLE 2 
MCMC Joint Linear Fits 



Sample 


A 


a 


Aft 


PAa.Ab 


X 2 /dof 


P X 2 


All galaxies 


-0.75 


+0.19 
-0.15 


30.4 


+6.9 
-7.2 


>99.9% 


7.3/4 


87.9% 


All galaxies (optimistic) 


-0.56 


+0.10 
-0.13 


23.5 


+5.1 
-4.5 


>99.9% 


18.9/5 


99.8% 


Norman et al. (2004) 


-0.38 


+0.09 
-0.11 


17.0 


+3.9 
-4.7 


>99.9% 


14.0/7 


94.9% 


Early-type galaxies 


-0.58 


+ 0.25 
-0.21 


21.6 


+ 11.4 
-7.7 


99.4% 


4.7/4 


68.4% 


Late-type galaxies 


-1.03 


+0.33 
-0.23 


39.1 


+ 12.8 
-9.9 


>99.9% 


3.5/3 


67.4% 


Early/late-type galaxies, low-z 


-0.19 


+0.15 
-0.14 


7 9 +6.0 
'• Z -5.6 


89.8% 


6.6/4 


84.4% 


Early /late- type galaxies, hi-z 


-0.59 


+0.32 
-0.34 


22.8 


+ 15.4 
-11.7 


96.9% 


1.6/3 


33.2% 



Note. — Best-fitting change in slope (Aa) and intercept (Aft) from joint fits to the low and high-z samples. The "Early /late- type" 
samples refer to comparing the early and late- type galaxy sample (i.e., see Figure [T2t . PAa Ab refers to the probability that Aa < and 
Aft > 0. 

P x 2 gives the \ 2 probability at which the model fit can be rejected (note that \ 2 is computed excluding upper limits). 
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TABLE 3 

Log-Normal Fits to the 60 [im Luminosity Function 



Parameter Saunders et al. (1990) Takeuchi et al. (2001) MCMC (This Work) 



a 


1.09 ± 0.12 


1.23 ±0.04 


1 04+ - 07 


a 


0.724 ±0.031 


0.724 ±0.010 


n 7^1 +0.014 
u -' oi -0.015 


<f> h 3 Mpc~ 3 


0.026 ± 0.008 


0.026 ± 0.003 


026+ 003 
U.UZD_ 003 


L* h~ 2 L Q 


10 8.47±0.23 


(4.34 ± 0.86) x 10 8 


O OQ+0.12 

10 -0.15 



Note. — The data used in our MCMC fitting of the 60 /jm LF 
was taken from Saunders et al. (1990), with several upper-limit 
data points excluded. 
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Fig. 14. — Log-normal (solid curves) and Schechter function (dashed curves) fits to the full (early + late- type) galaxy sample. The left 
panel shows the z < 0.5 sample and the right panel shows the 0.5 < z < 1.2 sample. 




Lt (ergs s l ) Lt (ergs s" 1 ) 

Fig. 15. — Log-normal (solid curves) and Schechter function (dashed curves) fits to the optimistic galaxy sample. The left panel shows 
the z < 0.5 sample and the right panel shows the 0.5 < z < 1.2 sample. 
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Fig. 16. — Log-normal (solid curves) and Schechter function (dashed curves) fits to the N04 galaxy sample. The left panel shows the 
z < 0.5 sample and the right panel shows the 0.5 < z < 1.2 sample. 




Lt (ergs s ') (ergs s ') 

Fig. 17. — Log-normal (solid curves) and Schechter function (dashed curves) fits to the early-type galaxy sample. The left panel shows 
the z < 0.5 sample and the right panel shows the 0.5 < z < 1.2 sample. 
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(ergs s ') L,, (ergs s ') 

Fig. 18. — Log-normal (solid curves) and Schechter function (dashed curves) fits to the late-type galaxy sample. The left panel shows 
the 2 < 0.5 sample and the right panel shows the 0.5 < z < 1.2 sample. 




Fig. 19. — Marginalized posterior probabilities for the fit parameters log0*, logL*, a and a for the log-normal fit to the z < 0.5 early-type 
galaxy XLF. The solid (black) lines show the posterior probability, the dotted (green) line show the prior, the dashed (red) lines shows 
a Gaussian distribution with the same mean and standard deviation as the posterior, and the dot-dashed line shows the 68% confidence 
interval. 
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log L' log L~ 

Fig. 20. — The z < 0.5 early-type (left) and late-type (right) galaxy MCMC draws for logci* and logL* in the log-normal fit. The 
solid lines show the 68% confidence intervals determined from the marginalized (and hence one-dimensional) posterior probabilities for the 
parameters. 



TABLE 4 
MCMC Log-Normal Fits 



Sample 


logcA* 


logL* 


a 


<7 


logp 


X 2 /dof 


P X 2 


All galaxies, low-z 


9 90+O.I8 
Z - ZO -0.39 


3Q S4+ - 57 
oy.o^t_ 43 


1 4 6 +0-20 1 


71 +0.08 t 
u ' ,1 -0.07 


37 94+0-07 


8.1 


N/A 


All galaxies (optimistic), low-z 


-2 46+ ' 31 


40 17+ - 44 
W.ll _Q.49 


1 46 +0-l7t 
1 '*°-0.19 


7 o+0.08 t 
u - '°-0.07 




1.8/1 


79.9% 


Norman et al. (2004), low-z 


9 SS+ ' 34 


40 40+0- 51 


i = 7+0.15 


74+0-08 1 
u - '^-0.07 


37.96t°;^ 


0.5/2 


21.7% 


Early-type galaxies, low-z 


o r.o+0.36 

o - vo -0.30 


40 13+ 48 


1 42+°- 14 t 
1,4 -0.19 


74+0-06 1 
u - ^-0.09 


37 62 +013 
°'- DZ -0.16 


3.0 


N/A 


Late-type galaxies, low-z 


-2 41+ - 34 


39.6lt°;£ 


1 46+0-35 1 


Q 7 o+0.07 t 


37.63±°;Jg 


1.6 


N/A 


All galaxies, hi-z 


-2 14+ - 48 
z - x -0.48 


39.85j:°;«g 


1 36 +0 45 t 


75+0-07-r 
u - '°-0.08 


38.27+g- 2 7 


0.9 


N/A 


All galaxies (optimistic), hi-z 


-2 34+ - 27 
z -° -0.31 


40 32+0- 52 
iu.oz_ Q 4g 


-, qo+0.37f 
1 " 38 -0.34 


Q 71 +0.08 t 
u - (± -0.07 


38 44+0-' 14 


1.2 


N/A 


Norman et al. (2004), hi-z 


-3 4+ 10 

° — 2 2 


4i.2+; : j 


O 09+ 016 t 
z.U3_ 27 


7 o+0.08 t 
u - ' J -0.07 


qo qq-I-0.56 

oo.zz_ 09 


6.0/1 


94.1% 


Early-type galaxies, hi-z 


o in + 0.54 




1 42+°- 2 °t 
1,4 -0.20 


Q 7 o+0.08 t 


37.83l°; 1 1 6 3 


1.6 


N/A 


Late- type galaxies, hi-z 


-1 76+ ' 81 


3Q 43+0-57 


, on+0.40 1 
1 - OZ -0.37 


72 +0.08 t 
u - ,z -0.06 


38.24i»;i| 


0.1 


N/A 



Note. — Best-fitting parameters from fitting a Log-Normal function to the XLFs. 
t Parameter is tightly constrained by prior 
* Parameter is moderately constrained by prior 
Luminosities are in ergs s — 1 in the 0.5-2.0 keV bandpass, 
p is in ergs s _1 Mpc -3 in the 0.5-2.0 keV bandpass. 

P x 2 gives the x 2 probability at which the model fit can be rejected (note that x 2 is computed excluding upper limits). 



TABLE 5 

MCMC SCHECHTER FITS 



Sample 


log</>* 


logL* 




Q 




logp 


X 2 /dof 


P X 2 


All galaxies, low-z 


q Qiy-|-0.28 

— o.z;_ 26 


40.931"; 2 * 


-1 


4^+0-16 
4o -0.12 


37 


qn+0.08 
MU -0.08 


2.9/1 


87.9% 


All galaxies (optimistic), low-z 


_3 45+0-24 
°-*°-0.27 


4l.28l°; 2 | 


-1 


ou -0.11 


38 


no+0.08 
08-0.07 


0.4/2 


16.9% 


Norman et al. (2004), low-z 


4 1S+ - 27 


41.74±g-.H 


-1 


66+ - 06 
— 0.09 


37 


qq+009 
yM -0.08 


2.1/3 


45.0% 


Early- type galaxies, low-z 


o =7+0.35 
c> - ' -0.36 


40.94l°;i 


-1 


,1 +0.24 
OA -0.22 


37 


55+0.13 


3.3/1 


89.7% 


Late-type galaxies, low-z 


o t-i+0.41 
°- ol --0.B7 


40.83lg;« 


-1 


KO+0.18 
OZ -0.25 


37 


fi o+0.11 


0.9/1 


62.9% 


All galaxies, hi-z 


_ 4 14+0.59 


41 71 + - 54 


-2 


04+0-26 1 
U4 -0.35 


38 


9 7+0.87 
z ' -0.21 


1.5/1 


75.7% 


All galaxies (optimistic), hi-z 


_3 54+0.25 


41.67lg' 2 ° 


-1 


R-i+0.23 
Di -0.25 


38 


49+0.20 
4Z -0.13 


0.3/1 


43.5% 


Norman et al. (2004), hi-z 


4 no+0.47 


42.05l°; 6 2 5 6 


-2 


no+0.28 
uo -0.12 


38 


97+0.41 
zl -0.14 


19.8/2 


99.5% 


Early-type galaxies, hi-z 


_4 31+O.53 


7 +0.74 
4i - ' '-0.32 


-1 


O4+0.41 i 
° 4 -0.32 


37 


77+0.42 
' ' -0.20 


4.2/1 


92.2% 


Late- type galaxies, hi-z 


_3 41+0.29 
J - 41 -0.37 


41 14+ ' 22 
^ 1 - i4 -0.17 


-1 


55+0.501: 
oo -0.38 


37 


on + 0.40 

o,:> -0.20 


0.9 


N/A 



Note. — Best-fitting parameters from fitting a Schechter function to the XLFs. 
t Parameter is moderately constrained by prior 
Luminosities are in ergs s — 1 in the 0.5-2.0 keV bandpass, 
p is in ergs s _1 Mpc -3 in the 0.5-2.0 keV bandpass. 

p x 2 gives the x 2 probability at which the model fit can be rejected (note that x 2 is computed excluding upper limits). 
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log/) (orgs s 1 Mpc" :: ) \ogi> (or^a s -1 Mpc 

Fig. 21. — The posterior probability distributions for the luminosity density (p) derived from log-normal fits to the early-type (left panels) 
and late-type (right panels) galaxy samples. The z < 0.5 results are shown in the top row and the 0.5 < z < 1.2 results are shown in the 
bottom row. The lines are as shown in Figure IT9l (there is no prior for p since it is derived from other fit parameters). 
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log/* (crga s 1 Mpc" :: ) k>gt> (orgs a -1 Mpc ,! ) 



Fig. 22. — Marginalized posterior probabilities for the fit parameters log</>*, logL*, and a for the Schechter function fit to the z < 0.5 
early-type galaxy XLF . Also shown is the probability distribution for the luminosity density p derived from the fit parameters. Lines are 
as shown in Figure Il9l 




Lx (ergs S ') L x (ergs s ') 

Fig. 23. — Joint fits to the low-z and high-z XLFs for the full (early + late-type) galaxy sample with the log-normal (left) and Schechter 
functions (right). In these fits the function parameters are tied between the two XLFs, and the offset in log L* is introduced as an additional 
fit parameter (i.e., pure luminosity evolution is assumed). The low-z XLF points are marked with (black) squares and the high-z XLF 
points are marked with (blue) triangles. The dashed (black) lines and dot-dashed (blue) lines show the fit to the low-z and high-z XLFs. 
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4. DISCUSSION 

We used the GOODS survey to derive X-ray luminos- 
ity functions for sources segregated by optical/NIR spec- 
tral type. We split the sources into low (z < 0.5) and 
high (0.5 < z < 1.2) redshift samples in order to inves- 
tigate evolution. We also explored an "optimistic" sam- 
ple (when our galaxy selection criterion is relaxed) and 
the N04 XLFs for comparison. We implemented MCMC 
techniques for linear, log-normal and Schechter function 
fits to the binned XLFs. This has given us a reliable sta- 
tistical assessment of the XLFs and any evolution. While 
in general either a log-normal or Schechter function could 
fit a given XLF well, this is due in part to the relatively 
sparse sampling in the XLFs presented here. Better data 
would be required to constrain the shapes of the XLFs. 
A consequence of this is that the faint-end slopes of the 
XLFs are somewhat uncertain, with the log-normal and 
Schechter function fits giving divergent predictions for 
the numbers of galaxies expected in deeper exposures. 

4.1. Comparison with Local XLFs for Different Spectral 

Types 

We showed that the difference in the early and late- 
type galaxy XLFs was only significant at the ~ 90% 
and - 97% level for the z~ 0.25 and z ~ 0.75 XLFs. 
This is suggestive that there is a difference between 
the spectral type-selected XLFs at each redshift but 
clearly better data will be required to strengthen this 
result. Early and late - type X LFs were also derived in 
iGeorgantopoulos et al.l (|2005l hereafter G05) for a sam- 
ple of galaxies at z< 0.2. In Figure [29] we plot our low-z 
early and late- type normal/starburst sample XLFs with 
the G05 XLFs also shown. There is good overall agree- 
ment between the two sets of XLFs, with the exception 
of the G05 early-type point at Lx ~ 10 39 ergs s _1 be- 
ing marginally higher than our corresponding point. The 
mean redshift of our low-z galaxy XLFs is ~ 0.3 for both 
the early-type and late-type samples, tentatively imply- 
ing that there has been little or no evolution between z 
- 0.1 and z - 0.3. However, - 40% of the G05 galaxy 
sample was comprised of CDF sources, most of which 
probably are also in our normal/starburst galaxy sam- 
ple. Having overlapping sources between our GOODS 
and the G05 samples would obviously dilute any differ- 
ence between the XLFs. 

4.2. Evolution 

In all of the samples discussed here the low-z and high- 
z XLFs differ at a confidence of > 99%, showing there 
is statistically-significant evolution with redshift. When 
log-normal and Schechter functions are fit to the indi- 
vidual XLFs, not surprisingly the resulting <j>* and L* 
values differ for a given sample, however the implied lu- 
minosity densities p do not in general differ significantly. 
In other words, our computed values of p are insensi- 
tive to whether the log-normal or Schechter functions 
are used, giving additional confidence in the resulting 
values. We also fit the low and high redshift XLFs simul- 
taneously with the log- normal and Schechter functions, 
with the main function parameters (i.e., Iog0*, logL* 
and the slopes) tied between the two XLFs but assum- 
ing pure luminosity evolution by fitting also for A logL*. 
The results for A log L* also do not depend strongly on 



whether the log-normal or Schechter functions are used. 
The corresponding A log p values were nearly identical to 
the A logL* values, which might be expected since logp 
depends linearly on logL* when all other parameters are 
fixed (and the differences observed between A log L* and 
A log p show the impact of jointly varying log <fi* and the 
slopes). 

The early and late-type galaxy AlogL* values were 
0.23 and 0.34 (from the log-normal fits), showing stronger 
evolution for late-type galaxies. This corresponds to 
early and late- type galaxies being factors of ~ 1.7 and 
~ 2.2 brighter, respectively, between z ~ 0.25 and 
z ~ 0.75. It is not clear whether the evolution ob- 
served in early-type galaxies is due to pa ssive evolution of 
the low-mass X-ray bi nary population (Ghosh & White 
120011 : iPtak et alll2001f ) or is due to enhanced star forma- 
tion as expected in the case of late-type galaxies. How- 
ever we note that our galaxy type selection is based on 
galaxy SED type rather than morphological type, and 
therefore these galaxies have red colors. This suggests 
that they are not actively star forming. The lack of 
any redshift dependence in the Fy/NIR flux ratio (see 
Appendix C) may be constraining to LMXRB evolution 
models sinc e LMXRB do not "turn on" instantenously as 
discussed in lGhosh fc White] (|2001h . however the scatter 
here is large. 

Pure luminosity evolution is often expressed as L*(z) = 
(1 + z) p L*(z — 0). With that parametrization and 
A log L* being measured between two redshifts Z\ and z 2 , 
P = i ,u A n L " f u s or p ~ 6.8AlogL* for z x ~ 0.25 

r log(l+z 2 )-log(l+zi) r 

and Z2 ~ 0.75. The AlogL* values then correspond 
to p cariy = X.h7±\Z and Mate = 2.33lg;g. The late- 
type galaxy evolution is consi stent with the FIR evolu- 
tion of p ~ 3. Also note that iGcorgakak is et aTl (2007) 
similarly found p ~ 2.4 for star-forming galaxies from 
the GOODS-N, using methods somewhat independent of 
those discussed here (although in both studies low L x , 
X-ray hardness, and X-ray/optical flux ratio were among 
the selection criteria). The full sample AL* was 0.29, not 
surprisingly intermediate to the early-type and late-type 
galaxy XLF values. The optimistic sample resulted in a 
A log L* value of 0.35, basically the same as the late-type 
galaxy value. However this may be somewhat coinciden- 
tal since the optimistic sample is most likely also intro- 
ducing low-luminosity AGN. However, any AGN activity 
is probably not dominating the near-IR-optical and the 
X-ray bandpass since otherwise the X-ray/optical ratios 
and/or the X-ray hardnesses would have resulted in an 
AGN classification (see Appendix A). Nevertheless the 
results of our analysis of the optimistic sample probably 
represents a reasonable limit to the maximum amount of 
evolution expected for the soft X-ray emission from nor- 
mal/starburst galaxies between z ~ 0.25 and z ~ 0.75. 
We also note that the luminosity densities inferred for 
the full sample XLFs are very similar to the luminosity 
densities inferred for the N04 XLFs, while the optimistic 
sample XLFs resulting in luminosity densities ~ 0.2 dex 
higher. Better data (either from deeper Chandra expo- 
sures or future X-ray missions) would of course result 
in smaller errors on the X-ray properties of the sources 
which in turn improve the classification probabilities, as 
well as increasing the number of sources populating the 
low- luminosity end of the XLFs. 
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L x (ergs s ') L x (ergs S ') 

Fig. 26. — As in Figure [23] for the early-type galaxy sample. 




L x (ergs s ') L x (ergs S ') 

Fig. 27. — As in Figure [23] for the late-type galaxy sample. 



TABLE 6 
MCMC Joint Log-Normal Fits 



Sample 


log0* 


logL* 


Q: 




A log L* 


A log p 


X 2 /dof 


P X 2 


All galaxies 

All galaxies (optimistic) 
Norman et al. (2004) 
Early-type galaxies 
Late-type galaxies 


9 9 n+0.16 

z - ZJ -0.29 
_ 2 47 +0.29 
z - 4 ' -0.13 

-2 Q6+ ' 29 
^.MD_ 24 

9 Q8+0 31 

-2.y»_ 2g 

9 no+0.18 
z - z °-0.29 


39.74+g;« 
40.14l°; 3 4 ° 
40-52±g:fl 
4O.08±g;|l 
39.43+. q 1 ;™ 


1-^J_0.30 
1-^_0.19 
1 60+ 010 

i- DU -o.io 

1 49+0.21J 
l- 4z -0.27 
-, 2? +0.34 t 
l- z '-0.25 


u -' z -0.07 

U.Oy_o.06 
n ?4 +0.08t 

n 7 n+0.07 t 

06g +0.08t 


n 98 +0.09 

n or+O.OS 
u "-"°-0.08 
25+ 08 

0.23i°;it 

o.34l°;J? 


0.29+H? 

n ocr + 0.08 
U.JO_ .08 

27+ ' 07 
u - z '-0.09 

0.23± ;H 
0.33±°£ 


8.2/3 
2.6/4 
2.4/6 
5.0/3 
3.4/2 


95.7% 
38.0% 
12.2% 
82.7% 
81.4% 



Note. — Best-fitting parameters from fitting a Log-Normal function jointly to the low and high-z XLFs, allowing only logL* to vary 
(i.e., assuming PLE). 

t Parameter is tightly constrained by prior 
t Parameter is moderately constrained by prior 
Luminosities are in ergs s _1 in the 0.5-2.0 keV bandpass. 
p is in ergs s — 1 Mpc — 3 in the 0.5-2.0 keV bandpass. 

p x 2 gives the x 2 probability at which the model fit can be rejected (note that x 2 is computed excluding upper limits). 
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TABLE 7 
MCMC Joint Schechter Fits 



Sample 


log <f>* 


logL* 


a 


A log L* 


Alogp 


X 2 /dof 


P X 2 


All galaxies 

All galaxies (optimistic) 
Norman et al. (2004) 
Early-type galaxies 
Late-type galaxies 


o 7n +0.19 
— 0.22 

Q ^0+°' 12 

°- ou -0.17 

-4 23+ - 17 
-4 26+ ' 30 

q £?o + 0.28 
°' DZ -0.24 


41 24+ ' 19 

41.32+g;" 

41.75ig;j; 
41 47+0.34 
_o.23 

40.91+°;?° 


-1 fis+°- 10 

J " DO -0.09 
1 cn+0.07 

^ -0.09 

-J.-O8_n.06 
, co+0.14 
-J-OO-0.13 
_1 KQ+0.16 
i - oy -0.13 


o 2q+°' 11 

u - zy -0.09 
S5+ 07 

24+ ' 07 

26+° 17 

S5+ 011 


0-34l ;° 8 9 

p, 9 O+0.08 

u - zci -0.08 

24+ ' 19 
u - z *-0.18 

S5+ - 11 
u.oo_ 14 


6.3/4 
1.0/5 
4.1/7 
5.0/4 
1.6/3 


82.0% 

4.1% 

23.4% 

70.8% 

35.0% 



Note. — Best-fitting parameters from fitting a Schechter function jointly to the low and high-z XLFs, allowing only log V to vary 
(i.e., assuming PLE). 

Luminosities are in ergs s - 1 in the 0.5-2.0 keV bandpass. 
p is in ergs s _1 Mpc - 3 in the 0.5-2.0 keV bandpass. 

P x 2 gives the \ 2 probability at which the model fit can be rejected (note that \ 2 is computed excluding upper limits). 
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Fig. 28. — The posterior probability distributions for AlogL* for the early-type (left) and late-type (right) galaxy samples based on log- 
normal fits (very similar results were obtained from the Schechter function fits). The solid solid (black) lines show the posterior probability, 
the dashed (red) lines shows a Gaussian distribution with the same mean and standard deviation as the posterior, the dotted (green) line 
show the prior, and the dot-dashed line shows the 68% confidence interval. 
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Fig. 29. — The low-z spectral type and galaxy selected XLFs 
from this work shown along with the XLFs from G05. 

5. SUMMARY 

We have computed XLFs for normal/starburst galaxies 
in the GOODS for sources with X-ray counterparts, and 
fit the XLFs with linear models using "traditional" tech- 
niques as well as Markov-Chain Monte Carlo techniques. 
From the photometric redshift fitting procedure we clas- 
sified 40 galaxies as early-type and 46 galaxies as late- 
type based on their SEDs. The early- type galaxy XLFs 
tend to be slightly flatter than those for late-type galax- 
ies, although from the MCMC analysis the significance 
of this result is only at the 1 — 2a level. The early and 
late-type galaxy sample XLFs at 0.25 are consistent 
with the low-redshift early and late-type XLFs of Geor- 
gantopoulos et al. (2005). We used the MCMC approach 
to also fit the XLFs with log-normal and Schechter func- 
tions. The XLFs discussed here all show significant evo- 
lution between z ~ 0.25 and z ~ 0.75. We jointly 
fit the low and high-redshift XLFs assuming pure lu- 
minosity evolution, allowing only L* to vary between 



the XLFs, which resulted in evolution of (1 + z) 16 and 
(l + z) 2 - 3 for the early-type and late- type galaxy samples. 
The late-type galaxy evolution derived here is consistent 
with the star-forming ga laxy X-ray evolution given in 
Gcorgaka kis et ail (|2007| ) . Including sources with am- 
biguous classification results in an "optimistic" galaxy 
sample with a total galaxy X-ray evolution of (1 + z) 2A , 
essentially the same value as derived for late-type galax- 
ies. The optimistic sample XLF evolution suggests that 
the maximum amount of evolution in the X-ray emis- 
sion of normal/starburst galaxies at these redshifts is 
(l + z) 2 - 4±0 - 5 . 

The Bayesian fitting approach here could be expanded 
to include additional uncertainties that might impact 
this analysis , such as the redshift errors (see, e.g., 
iDahlen et alj [2005) , the uncertainties in the complete- 
ness correction and the uncertainties the X-ray and op- 
tical fluxes of the sources. A larger impact on our re- 
sults would likely result from including radio and FIR 
data from the GOODS fields, which we will explore 
in future work. This will improve both the SED fit- 
ting (and hence the galaxy type determination) and give 
an independent star-for mation rate estimation (see also 
iGeorgakakis et al1l2007f ). It may also be possible to si- 
multaneously fit for the multi-variate luminosity func- 
tions and the individual galaxy types, activity types and 
redshifts (with any spectroscopic redshifts used as tight 
priors), at least in an iterative fashion (i.e., where the 
current luminosity function estimates guide the galaxy 
type and photometric redshift probabilities). Finally, 
advanced Bayesian mod el selection techiniques, as dis- 
cussed in iTrottal {2007) and references therein, can be 
applied here to guide the parameters of future observa- 
tions by predicting the ability of future data to prefer a 
given model and arrive at a given set of constraints. 
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APPENDIX 

GOODS X-RAY SAMPLE PROPERTIES AND GALAXY SELECTION 

Here we discuss our methodology for classifying the sources, and we also discuss the statistical properties of the 
sample. Note that the relevant point of our classification is not whether the entire SED is dominated by star-formation 
or an AGN, but which of these is dominating the soft X-ray band. 

Bayesian Selection 

N04 selected normal/starburst galaxies from the full CDF-N and CDF-S samples using a Bayesian classification pro- 
cedure, where priors were constructed from the a set of galaxies with well-determined optical types, normal/starburst 
galaxy, type-1 AGN and type-2 AGN (hereafter galaxy, AGN1 and AGN2). The product of the prior distributions for 
a class and the likelihood for the observed parameters for a given source gave the probability that the source was drawn 
from that class. Here we follow the same procedure with several improvements. First, sources with small differences 
between the probabilities in each class should have been considered uncertain since these conditions occur when the 
separation of the sources parameter values and the parent distribution means are small relative to the parameter 
errors. Here we use the Bayesian "odds ratio" , or the ratio of posterior probabilities for the classes being compared. 
In Bayesian model testing a model is considered to be favored only when the odds ratio exceeds at least 3 (while 
odds ratios greater than 10 are preferred). Second, the parameter likelihoods were handled somewhat simplistically, 
with Gaussian errors assumed on each parameter and a constant error in \ogLx and logL^ — logL opt of 0.25. The 
Gaussian assumption is not correct for hardness ratio errors and logLx errors when the number of counts detected 
for the source is s mall. However, it turns out that the Gaussian approximation to hardness ratio errors is conservative 
(|Park et al.ll2006fh Here we use the larger of the asymmetric errors on count rate given in Alexander et al. (2003). 
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We assume that AFx /Fx = AC/C where AFx is the error on X-ray flux and AC is the error on X-ray count rate 
C given in Alexander et al. We then take the error on logLx to be given by L J\^ 10 ) where ALx = Lx A p* ■ This 
is valid only when AFx I Fx <C 1, however this is the case for the majority of the X-ray sources. Finally, N04 did 
not account for k-corrections in the optical data when computing X-ray /optical flux ratios (k-corrections to the X-ray 
data were not necessary since an energy index of 1.0 was assumed for every source). As discussed below, k-corrections 
are now included in the computation of the priors and in the galaxy classification. 

Priors 

Szokoly et al. (2004) reported classifications based on the optical spectra alone, in the classes "ABS" (no or only 
absorption lines are present in the spectrum), "LEX" (a low-ionization emission line spectrum), "HEX" (a high- 
ionization emission line spectrum), and "BLAGN" (broad emission lines are found). All of the BLAGN sources should 
correspond to AGN1 sources, by definition, however a broad line AGN may be present in the other classes where 
there was not sufficient signal to detect a broad-line component. The HEX and LEX classes should be dominated 
by type-2 AGN and star-forming galaxies, respectively. However the LEX cl assification includes some AGN where 
low signal-to-noise or dilution of the nuclear spectrum due to aperture effects (|Moran et alJl2002T ) has precluded the 
detection of high-excitation emission lines. For several LEX sources the statistics are sufficient for an AGN component 
to be identified from line-ratio diagnostics. We therefore derive priors using only ABS sources as galaxies, HEX sources 
and LEX sources with AGN line ratios as AGN2, and BLAGN sources as AGN1. We required that the X-ray sources 
have a corresponding entry in the Alexander et al. catalog since that catalog is used for the X-ray properties of the 
GOODS sample. 

We initial only selected sources with z < = 1.2 however this resulted in only 8 AGN2 sources. Since AGN2 are known 
to have flat X-ray spectra, and hence a minimal k-correction, we relaxed the redshift constraint to z <= 2 for that 
class. The final tally was then 11 AGN2 sources, 11 AGN1 sources and 15 (normal/starburst) galaxy sources. The 
mean offsets between the R and K band magnitudes reported in Szokoly et al. and the R and K s bands used in the 
GOODS survey were computed in order to adopt the priors based on the Szokoly et al. source for use with GOODS 
data. This was done regardless of spectral type and we found offsets of 0.22 magnitudes in R and 1.9 magnitudes in 
K, with a standard deviations of 0.16 and 0.21 magnitudes. 

As shown in Figure [3] k-corrections to the X-ray/R and X-ray/K band flux ratios can be significant. Since only K 
and R band magnitudes are listed in Szokoly et al., we cannot apply the same interpolation procedure for k-corrections 
as was done with GOODS data. However, we can use the R-K color to estimate the spectral type of the source, and 
apply a mean k-correction based on that type. We show in Figure [30] the R-K s colors for the full sample, with early- 
type, late-type and irregular/starburst galaxies marked separately. While there is some overlap, we manually selected 
linear functions to delineate the spectral type as a function of redshift, as shown in the figure. The ear ly/late- type 
galaxy separation is given by R-K s = 1.0 + 2.2z and the late- type/star burst galaxy separation is given by R-K s = 
0.2+2. Oz. 

Having established a crude spectral type color selection, we proceeded by plotting the R and K s k-corrections 
as a function of redshift for each spectral type, shown in Figure 1311 For each spectral type we determined the 
k-correction/redshift correlation. Again, while the dispersion in k-correction is large, this resulted in at least an 
approximate k-correction that was applied to the sources when computing the priors. 

We list in Table [8] the mean and standard deviation for the parameters Lx, HR (X-ray hardness), X-ray/R-band 
flux ratio and X-ray/K s -band flux ratio computed using the prior samples. The priors were then modeled as Gaussians 
with these values for the Gaussian mean and standard deviation, as discussed in N04 The full sample is listed in Tables 
[TT1and[C| where we list the final classification along with the Bayesian odds ratios for the models. For example, if the 
preferred class for a source is "galaxy", then the odds ratios are given f or ga laxy vs. AGN1 and galaxy vs. AGN2. 
For comparison, 30/43 of the OBXF sources from iHornschemeier et all (|2003l ) are in our GOODS-N sample, and all 
but one were classified as (normal/starburst) galaxies. The remaining source (XID=458) had "galaxy" as the highest 
probability, but only at a factor of 2.1 higher than the AGN2 probability and hence is only in the "optimistic" sample. 
A spectral type was given for 28 of the 30 OBXF sources with matches, which is also listed in Table [TT1 Overall there 
was good agreement between our SED types and the spectral types, with 7/9 of the absorption-line galaxies having 
an early- type SED classification and 15/17 emission-line galaxies having a late- type SED classification (the remaining 
2 sources were "composite" galaxies having both absorption and emission lines). Similarly we also list the redshift 
found in Georgakakis et al. (2007), if present, for the GOODS-N sources, where it can be seen there is good overall 
agreement between the redshifts. 

STATISTICAL PROPERTIES OF THE SAMPLES 

In Figure [32] we show histograms of the photon indices, binned separately by spectral type. The photon index 
distributions are similar, peaking at T ~ 1.2 — 1.4. However note this is due in part to the adoption of T = 1.4 
in Alexander et al. for sources with low signal-to- noise. The number of irregular/starburst sources is very low and 
is likely to be similar to the rest of the sample. We list in Table [5] the number of sources in each sample (i.e., 
divided by field, host galaxy type and bandpass) along with the mean, standard deviation, minimum and maximum 
(k-corrected) magnitude. We also give the similar statistics for log Fx and log Lx ■ Note that K s magnitudes were not 
always available. An immediate conclusion is that the AGN contribution is most significant for the irregular/starburst 
samples since the mean X-ray luminosities are ~ 10 42 ergs -1 , the luminosity at which AGN emission starts to dominate 
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TABLE 8 
Bayesian Prior Parameters 



Class 


Parameter 


Mean 


17 


Normal/Star burst Galaxies 


HR 


-0.19 


0.46 




log L x 


40.6 


0.7 




log F x /F opt 


-3.2 


0.7 




log F x /F K 


-3.4 


0.7 


AGN2 


HR 


0.16 


0.37 




\ogL x 


41.1 


1.1 




log F x /F opt 


-2.2 


1.0 




log F x /F K 


-2.7 


0.7 


AGN1 


HR 


-0.51 


0.05 




log L x 


42.9 


0.4 




log F x /Fo P t 


-1.2 


0.4 




log F X /F K 


-1.4 


0.5 
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Fig. 31. — K-corrections plotted as a function of redshift for R-band (left) and K s -band (right) magnitudes. The quantity plotted is the 
difference between the observed and k-corrected magnitude. Early-type galaxies are plotted with diamonds , late-type galaxies are plotted 
with pluses, and irregular/starbursts are plotted with asterisks. 
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Fig. 32. — Distribution of photon indices for the CDF-N (left) and CDF-S (right) sources, binned by galaxy spectral type (black, solid 
line=early-type, red, dotted line=late-type, blue, dashed line=starburst/irregular). 



the X-ray band. The corresponding values after galaxy selection are listed in Table [TO] 
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TABLE 9 

X-ray and Optical Statistical Properties 



Field 


Band 


SED Type 


N 


Mean 


a 


Min. 


Max. 


CDFN 


B 


ear y 


57 


22.26 


1.07 


19.13 


24.27 


CDFN 


R 


ear y 


57 


21.29 


1.16 


18.03 


23.45 


CDFN 




early 


^7 


io nn 
lo.yu 




17 nn 


21 32 


CDFN 


J 


early 


57 


20.32 


1.02 


17.45 


22.32 


CDFN 


logFX 




57 


-16.16 


0.48 


-17.38 


-14. 64 


CDFN 


logLX 


early 


57 


40.96 


0.70 


38.91 


42.42 


CDFN 


B 


late 


73 


21.94 


1.16 


19.64 


24.66 


CDFN 


R 


late 


73 


21.20 


1.15 


18.59 


23.81 


CDFN 


Ks 


late 


73 


19.18 


0.85 


17.62 


21.69 


CDFN 


j 


late 


73 


20.41 


0.95 


18.09 


22.53 


CDFN 


logFX 


late 


73 


-16.03 


0.79 


-17.39 


-13.64 


CDFN 


logLX 


late 


73 


41.01 


1.03 


38.89 


43.33 


CDFN 


B 




g 


22.11 


1.48 


19.91 


24.54 


CDFN 


R 




g 


21.73 


1.38 


19.84 


23.88 


CDFN 


Ks 


ir 


g 


20.20 


1.55 


18.27 


22.47 


CDFN 


j 




g 


21.18 


1.39 


19.69 


23.13 


CDFN 


logFX 


! r 


g 


-15.27 


1.27 


-16.74 


-13.81 


CDFN 


lOgljvY 




Q 
O 


ZL1 7^ 


1 52 


79 




CDFS 


B 


ir i 

ear y 


36 


21.62 


1.48 


18.11 


25.95 


CDFS 


R 


ear y 


36 


20.54 


2.08 


16.25 


27.61 


CDFS 


Ks 


ear y 


29 


18.73 


1.22 


16.07 


21.59 


CDFS 


j 


ear y 


29 


19.86 


1.38 


16.27 


22.80 


V .. 1 I i o 


logr <-\ 


early 


0\J 


- lO.OU 


u.ou 


1 £8 
-1U.OO 


-14 51 


CDFS 


logLX 




36 


41.32 


0.65 


39.65 


42.77 


CDFS 


B 


late 


32 


22.03 


1.68 


17.34 


25.95 


CDFS 


R 


late 


32 


21.08 


1.70 


16.71 


24.97 


CDFS 


Ks 


late 


26 


19.08 


1.27 


15.94 


21.82 


CDFS 


J 


late 


26 


20.25 


1.50 


16.03 


22.39 


CDFS 


logFX 


late 


32 


-15.83 


0.63 


-16.74 


-14.39 


CDFS 


logLX 


late 


32 


41.31 


0.83 


39.91 


42.90 


CDFS 


B 


ir 


13 


21.45 


1.55 


19.24 


24.65 


CDFS 


R 


ir 


13 


21.00 


1.69 


18.00 


24.30 


CDFS 


Ks 


ir 


7 


20.76 


1.33 


19.58 


22.89 


CDFS 


J 


ir 


7 


21.45 


1.30 


19.98 


23.64 


CDFS 


logFX 


ir 


13 


-14.90 


1.13 


-16.27 


-13.36 


CDFS 


logLX 


ir 


13 


41.98 


1.78 


39.44 


44.03 
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TABLE 10 

X-ray and Optical Statistical Properties , Galaxies only 



Field 


Band 


SED Type 


N 


Mean 




Min. 


Max. 


CDFN 


B 


early 


27 


21.54 


0.91 


19.13 


22.91 


CDFN 


R 


early 


27 


20.50 


1.00 


18.03 


21.88 


OUr IN 


rvS 


early 


Z I 


lo.ol 


U.ol 


1 7 no 

i / .uy 


iy.00 


OUr IN 


T 

J 


early 


2 / 


19.60 


0.81 


1 1 .45 


20.75 


CDFN 


logFX 


early 


27 


-16.32 


0.31 


-16.93 


-15.60 


UUr IN 


logLA 


early 


2 / 


40.63 


0.66 


38.91 


A 1 A *7 
41.4 / 


CDFN 


B 


late 


36 


21.45 


1.01 


19.64 


23.65 


OLtr IN 


K 


late 


36 


20.73 


1.04 


18.59 


22.87 


OUr ^N 


AS 


late 


OD 


lo.yo 


U.DD 


It . to 


on c^o 


CDFN 


J 


late 


36 


20.02 


0.80 


18.09 


21.45 


ODr IN 


1 T7 V 

logr A 


late 


36 


-16.39 


0.32 


-17.12 


-15.72 


CDFN 


logLX 


late 


36 


40.41 


0.64 


38.89 


41.41 


CDFN 


logFX 


ir 


1 










UUr ^N 


logLA 


ir 


1 










ODr b 


T> 
D 


early 


13 


20.47 


1.01 


18.11 


21.52 


ODr b 


o 
K 


early 


13 


19.15 


1.38 


1 a ok 
lb. 25 


20.96 


OD1 b 


Ks 


early 


11 


17.63 


0.76 


16.07 


18.74 


OD1 b 


J 


early 


11 


18.58 


1.07 


16.27 


20.09 


OL)r b 


1 T7 V 


early 


13 


-15.97 


0.46 


-16.59 


-15.18 


OD1 b 


logLX 


early 


13 


40.92 


0.62 


39.65 


41.85 


CDFS 


B 


late 


10 


20.52 


1.61 


17.34 


21.83 


CDFS 


R 


late 


10 


19.70 


1.58 


16.71 


21.23 


CDFS 


Ks 


late 


9 


18.08 


1.21 


15.94 


19.54 


CDFS 


J 


late 


9 


18.92 


1.56 


16.03 


20.11 


CDFS 


logFX 


late 


10 


-15.94 


0.43 


-16.36 


-15.07 


CDFS 


logLX 


late 


10 


40.62 


0.49 


39.91 


41.43 


CDFS 


B 


ir 


2 


21.03 


0.79 


20.47 


21.59 


CDFS 


R 


ir 


2 


20.63 


0.74 


20.11 


21.15 


CDFS 


Ks 


ir 


2 


19.66 


0.11 


19.58 


19.74 


CDFS 


J 


ir 


2 


20.32 


0.48 


19.98 


20.66 


CDFS 


logFX 


ir 


2 


-16.24 


0.05 


-16.27 


-16.20 


CDFS 


logLX 


ir 


2 


40.35 


0.48 


40.01 


40.69 
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Fig. 33. — The X-ray/R-band flux ratio plotted as a function of X-ray luminosity, including (left) and not including (right) k-corrections. 
Early- type, late- type and irregular/starburst SED sources are marked with (black) diamonds, (red) pluses, and (blue) asterisks, respectively. 
The ellipses show the 1 and 2a intervals computed for the prior probability distributions for galaxies, AGN1 and AGN2 for logLx and 
log Fx I Fr- Note that these intervals correspond to standard deviations computed separately for these parameters (i.e., the widths of the 
ellipses give the 68% and 95% probability intervals for log Lx for each type) , which assumes no correlation between L x and Fx / Fr . 
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Fig. 34. — The X-ray/K s -band flux rati o plotted as a function of X-ray luminosity, including (left) and not including (right) k-corrections. 
Symbols and ellipses are as in Figure l33l 

X-RAY FLUX RATIOS 

Since X-ray/R-band and X-ray/K s -band flux ratios are used in selection criteria for this paper, we discuss here 
the potential impact of any luminosity or redshift dependence of these quantities. In Figures [33] and [M] we plot the 
X-ray/R-band and X-ray/K s -band flux ratios as a function of luminosity, before and after applying k-corrections. 
Figures [35l and [36l show the flux ratios plotted as a function of redshift. In the luminosity versus flux ratio plots we 
show the regions resulting from our Bayesian prior analysis. These ellipses show the 1 and 2a regions where the a 
value is based on the standard deviation of the parent populations. In other words, the width and height of the la 
regions are the standard deviations of the X-ray luminosity and given flux ratio for the corresponding class (galaxy, 
AGN1 or AGN2), and correspond to the 68% and 95% probability intervals for that parameter. However, this assumes 
that the X-ray luminosity and flux ratios are not correlated, and since this is not the case, these regions should not 
be interpreted as joint confidence regions. The intention here is to simply show the regions in the logLx — log Fx /Fr 
and logLx — log Fx /Fk b planes from which our normal/starburst samples are being selected. 

In the redshift versus flux ratio plots we show the la regions for the flux ratios, as in luminosity /flux ratio plots, 
except in this case the regions are simply marked with horizontal lines (i.e., since redshift is not a selection criterion). 
From these plots we conclude that any evolution in the flux ratios with redshift is dominated by the k-corrections 
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Fig. 35. — The X-ray/R-band flux ratio plotted as a function of redshift, including (left) and not including (right) k-corrcctions. Early- 
type, late-type and irregular /starburst SED sources are marked with (black) diamonds, (red) pluses, and (blue) asterisks, respectively. 
Horizontal lines are drawn showing the la probability intervals for the flux ratio for galaxies (green dashed lines), AGN1 (red dotted lines) 
and AGN2 (blue solid lines). 
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Fig. 36. — The X-ray/K s -band flux ratio plotted as a function of redshift, including (left) and not including (right) k-corrections. Early- 
type, late-type and irregular /star burst SED sources are marked with (black) diamonds, (red) pluses, and (blue) asterisks, respectively. 
Horizontal lines are drawn showing the la probability intervals for the flux ratio for galaxies (green dashed lines), AGNf (red dotted lines) 
and AGN2 (blue solid lines). 

and the scatter in the flux ratios. In other words, after k-correction, the flux ratios are consistent with no redshift 
dependence. This emphasizes the importance of k-correcting the data prior to utilizing flux ratios as a selection 
criterion. On the other hand there is clearly a luminosity dependence in the flux ratios. This is due in large par t 
to the increased prevalence of AGN activity at higher X-ray luminosities (e.g.. iFiore et alj|2003t iBarger et alj|2005j) . 
The correlation between X-ray luminosity and the flux ratios mitigates the benefit of including both luminosity and 
the ratios as selection criteria. However, in practice there are a significant number of sources in the luminosity range 
i()4i-42 ergs s _1 , consistent with either AGN or normal/star burst galaxies, where the X-ray/R-band and/or X-ray/K s - 
band flux ratio is within the la region for AGN2 galaxies. Thus, including the X-ray/optical and X-ray/K-band flux 
ratios as a selection criteria should improve the separation, and hence selection, of normal/starburst galaxies from 
type II AGN. 
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TABLE 11 
CDF-N X-ray Sample 



XID 


z 


Flag 




Class 


log(oddsl) 


log(odds2) 


log(Lx) 


log(F x /F R ) 


log(F x /F K ) 


Type 


OBXF ID 


ZG07 


48 


1.01 


4 


0.27 


agn 


38.53 


0.90 


42.16 


-1.97 


-2.99 


early 






55 


0.64 


\ 


0.12 








41.30 


-2.78 


-3.02 


late 






c^ft 


U. lo 


A 
4 


n on 
U.zU 


gal 


i/i ift 

14.10 


n £i 


on Q/i 

oy.y4 


O 70 
-Z. / Z 


o Qn 
-z.ou 


late 




n nfiat 
U.Uo 1 


57 


0.38 


2 


0.15 


gal 


12.78 


1.36 


40.87 


-3.21 


-3.33 


late 


hi (U.oo) 


0.33 a t 


60 


0.42 


4 


0.05 








40.54 


-2.42 


-2.52 


early 






62 


0.22 


4 


0.03 


gal 


18.12 


1.21 


39.58 


-4.56 


-4.44 


early 


A (0.09) 




bi 


0.64 


1 


0.14 


gal 


8.68 


0.98 


41.39 


-2.87 


-3.30 


late 


0.64 a 


1 Z 


n cia 


i 

1 


n 1 1 
U.ll 








A 1 ftQ 
41. DO 


O Qft 

-z.oD 


O 07 

-z.y / 


late 






78 
/O 


n it* 
U. / o 


i 

i 


n 1 q 
U. lo 


agnz 


i.oy 


ii/i 

1 . 14 


A 1 ftft 
41. DD 


1 70 
-1. 1 Z 


0/1 
-Z.Z4 


early 




n 7t^ a 
u. i o 


Q1 
ol 


n 

U.oo 


A 
4 


n 1 n 
U. 1U 


gal 


i n oo 
lu.yz 


n fi/i 
U.o4 


a n 7n 

4U. 1 U 


q in 
-o. 1Z 


OC 

-o.zo 


early 








n ftc 
U.Do 


1 

1 


n 1 a 

U. 14 








A 1 AA 
41.44 


O 

-Z.DO 


O Q7 
-Z.o / 


early 






£7 


n i /i 

U. 14 


1 

1 


n 1 o 

U. 1Z 


gal 


on ^7 

ZU.O l 


1 Q1 

l.ol 


on 7ft 

oy. (o 


-o.oy 


q on 
-o.zu 


early 






on 
yu 


11/1 


1 

1 


n 1 ft 
U. ID 


agn2 


1 Q 1 Q 

lo. lo 


n fi/i 
U.o4 


a o nft 
4Z.UO 


o no 
-z.uy 


-z.oo 


late 






yo 


U.zo 


1 

1 


k no 
o.Uo 


agn2 


1 on 
l.oU 


O 1 o 

z. Iz 


a o no 
4z.Uo 


i ?n 
-1. i u 


1 7c; 

-1. / 


late 






1U1 


n /i 
U.4o 


1 

1 


n ns 
U.Uo 


gal 


i o ftn 
Iz.DU 


1 1 o 

i.iy 


/in si 
4U.ol 


o. ns 
-o.Uo 


o Qn 
-o.oU 


early 


A (0.45) 


n a 

U.40 


103 


0.97 


2 


0.28 


agn2 


1.51 


0.66 


42.13 


-2.05 


-2.62 


late 




0.81 a t 


105 


0.33 


4 


0.03 


gal 


14.66 


0.96 


40.03 


-3.31 


-3.41 


early 






110 


1.01 


4 


0.34 


agn2 


0.89 


1.05 


42.26 


-1.98 


-2.54 


early 






111 


0.52 


1 


0.04 


gal 


13.33 


1.15 


40.64 


-3.17 


-3.57 


late 


E (0.52) 


0.52 a 


113 


0.84 


2 


1.59 


agnl 


4.98 


1.93 


42.74 


-1.29 


-1.90 


late 


0.84 a 


114 


0.53 


1 


0.04 


gal 


17.54 


0.65 


40.63 


-3.22 


-3.53 


early 






115 


0.68 


1 


4.67 


agnl 


5.48 


1.43 


42.98 


-1.21 


-1.38 


late 




0.68 a 


119 


0.47 


1 


0.08 


gal 


12.31 


1.14 


40.83 


-3.11 


-3.35 


early 


E (0.47) 


0.47 a 


120 


0.69 


1 


0.09 


gal 


7.17 


0.51 


41.27 


-2.78 


-2.75 


ir 




126 


0.77 


4 


0.03 


gal 


10.68 


0.93 


40.90 


-3.01 


-3.60 


late 




0.84 a 


132 


0.64 


1 


0.07 


gal 


11.36 


0.67 


41.06 


-3.01 


-3.26 


late 


E (0.65) 


0.64 a 


136 


0.47 


1 


0.05 


gal 


13.55 


1.16 


40.65 


-3.20 


-3.37 


early 


A (0.47) 


0.47 a 


138 


0.48 


1 


0.07 


gal 


10.81 


0.80 


40.80 


-2.92 


-3.05 


late 


E (0.48) 


0.48 a 


139 


0.93 


4 


0.76 


agnl 


2.87 


0.75 


42.52 


-1.60 


-2.20 


late 


1.01 a 


142 


0.75 


1 


0.29 


agn2 


7.05 


0.97 


41.87 


-2.02 


-2.60 


early 






150 


0.63 


4 


0.13 


agn2 


23.99 


1.38 


41.35 


-1.85 


-2.26 


early 






158 


1.01 


2 


0.17 


agn2 


20.44 


0.49 


41.96 


-2.16 


-3.33 


early 






160 


0.82 


4 


0.12 








41.57 


-2.73 


-2.91 


late 






166 


0.46 


1 


0.04 








40.50 


-2.78 


-2.63 


late 




0.46 a 


169 


0.31 


4 


0.10 


gal 


13.84 


1.10 


40.48 


-3.08 


-3.09 


late 


A (0.84) 




170 


0.63 


4 


0.20 


agn2 


11.44 


0.77 


41.53 


-2.13 


-2.62 


early 




177 


1.02 


1 


0.25 


agn2 


2.73 


0.66 


42.14 


-2.07 


-2.99 


late 




1.01 a 


180 


0.46 


1 


0.23 


gal 


9.37 


1.05 


41.25 


-2.97 


-3.35 


early 




0.46 a 


187 


0.94 


4 


0.04 


gal 


14.53 


0.49 


41.23 


-2.80 


-3.68 


late 






188 


1.15 


5 


0.07 


agn2 


4.05 


0.88 


41.71 


-1.93 


-2.47 


late 






189 


0.41 


1 


0.17 








41.01 


-2.57 


-2.71 


early 






194 


0.56 


1 


1.62 


agnl 


6.01 


2.94 


42.30 


-1.19 


-1.41 


late 




0.56 a 


197 


0.08 


1 


0.05 


gal 


28.79 


0.90 


38.90 


-3.81 


-3.43 


late 


E (0.08) 


0.08 a 


200 


0.97 


1 


0.04 








41.31 


-2.72 


-3.18 


late 




203 


1.14 


2 


0.03 


gal 


10.14 


0.98 


41.27 


-3.02 


-3.96 


late 






209 


0.51 


2 


0.23 








41.36 


-2.13 


-2.47 


late 




0.51 a 


91 n 
ziu 


n 7n 
u. / u 


/i 

4 


n oq 
u.uy 


gal 


0.40 


n 87 


A1 ^n 

41. OU 


71 
-z. / 1 


^ *3n 

-O.OU 


early 






211 


0.76 


4 


0.06 


gal 


8.76 


0.51 


41.23 


-2.82 


-3.10 


late 




0.85 a 


212 


0.94 


1 


0.13 






41.77 


-2.28 


-2.76 


early 






214 


1.04 


4 


0.03 


gal 


8.63 


0.75 


41.27 


-2.93 


-3.63 


early 






215 


1.01 


2 


0.03 








41.17 


-2.22 


-3.63 


late 






217 


0.54 


4 


0.04 


gal 


24.31 


0.79 


40.63 


-3.38 


-3.97 


early 






218 


0.09 


1 


0.10 


gal 


21.17 


0.65 


39.31 


-3.05 


-2.65 


late 


E (0.09) 




219 


0.70 


4 


0.03 


gal 


10.65 


0.85 


40.76 


-3.02 


-3.64 


early 




222 


0.77 


4 


0.59 


agn2 


3.17 


1.29 


42.21 


-1.88 


-2.60 


early 




0.86 a 



TABLE 11 — Continued 



XID 


z 


Flag 


Fx 


Class 


log(oddsl) 


log(odds2) 


log(Lx) 


log(F x /F R ) 


log(F x /F K ) 


Type 


OBXF ID 


ZG07 


227 


0.52 


4 


0.04 


gal 


13.17 


1.03 


40.57 


-3.14 


-3.34 


late 


E (0.56) 


0.56 a 


230 


1.01 


1 


0.08 


gal 


9.67 


1.11 


41.62 


-2.99 


-3.80 


early 




1.01 a 


234 


0.45 


1 


0.09 


gal 


10.56 


0.84 


40.85 


-2.89 


-2.97 


late 




0.45 a 




u.y / 


i 
i 


u.uo 


gal 


q Qn 
y.yu 


u.y i 


41.JO 


-O.UO 


-O.D 1 


late 




n Q7 a 
u.y / 


245 


0.32 


1 


0.02 


gal 


15.27 


0.81 


39.90 


-3.32 


-3.09 


late 


E CO 32~1 


0.32 a 


249 


0.47 


2 


0.24 








41.30 


-2.53 


-2.71 


early 




251 


0.14 


1 


0.07 


gal 


22.92 


1.31 


39.59 


-3.59 


-3.24 


late 


E (0.14) 


0.14 a 


256 


0.60 


1 


0.09 








41.10 


-2.65 


-2.72 


late 




257 


0.09 


1 


0.04 


gal 


30.11 


1.08 


38.92 


-4.09 


-3.93 


early 


A (0.09) 




258 


0.75 


1 


0.04 


gal 


9.97 


0.79 


40.96 


-2.83 


-3.34 


late 


0.75 a 


260 


0.47 


1 


0.06 


gal 


10.82 


0.86 


40.71 


-2.69 


-3.13 


late 




0.47 b 


262 


0.87 


4 


0.67 


agnl 


4.47 


1.93 


42.39 


-1.31 


-1.91 


early 




0.81 bt 


264 


0.32 


1 


0.05 








40.25 


-2.54 


-2.29 


late 




0.32 b 


265 


0.41 


1 


0.09 


p"al 


12.74 


1.30 


40. 7 4 


-3.03 


-3.35 


early 


C (0.41) 


0.41 b 


266 


1.08 


4 


0.04 








41.38 


-2.36 


-3.67 


early 




269 


0.36 


2 


0.05 








40.36 


-2.56 


-2.42 


late 




0.36 b 


274 


0.32 


1 


0.24 


gal 


10.66 


1.13 


40.91 


-3.01 


-3.13 


early 


A (0.32) 




278 


1.02 


2 


0.03 








41.24 


-2.49 


-3.16 


late 




279 


0.89 


2 


0.03 








41.07 


-2.46 


-2.92 


ir 






280 


0.96 


4 


0.04 








41.25 


-2.23 


-2.90 


early 






282 


0.20 


2 


0.06 


gal 


18.94 


1.10 


39.87 


-3.42 


-3.42 


late 


E (0.20) 


0.08 bt 


286 


0.95 


1 


0.41 


agn2 


0.71 


1.16 


42.28 


-1.92 


-2.61 


early 




0.95 b 


288 


0.71 


4 


0.11 


gal 


8.51 


0.88 


41.38 


-2.61 


-3.54 






0.79 b 


291 


0.52 


1 


0.04 


gal 


14.27 


1.20 


40.60 


-3.35 


-3.62 


early 






292 


0.50 


2 


0.02 








40.34 


-2.27 


-2.46 


ir 






294 


0.47 


2 


0.31 


agn2 


1.29 


1.50 


41.42 


-1.50 


-1.93 


late 




0.47 b 


295 


0.85 


1 


0.06 








41.31 


-2.36 


-2.68 


late 






296 


0.66 


1 


0.07 


gal 


11.67 


0.55 


41.12 


-2.89 


-3.30 


early 






300 


0.14 


2 


0.07 


gal 


25.94 


1.48 


39.54 


-3.91 


-3.75 


late 


E (0.14) 


0.14 b 


304 


0.68 


1 


3.22 


agnl 


5.87 


2.60 


42.81 


-1.30 


-1.75 


late 




305 


0.30 


1 


0.05 


gal 


17.69 


1.30 


40.15 


-3.43 


-3.47 


late 


E (0.30) 


0.30 b 


309 


1.14 


2 


0.68 


agnl 


4.25 


1.77 


42.70 


-1.37 


-2.28 


late 




1.14 b 


310 


0.76 


1 


0.03 


gal 


11.44 


1.02 


40.91 


-3.05 


-3.64 


late 




0.76 b 


311 


0.91 


1 


0.05 


gal 


7.15 


0.49 


41.33 


-2.56 


-3.18 


late 




0.91 b 


313 


0.80 


1 


0.07 


gal 


9.70 


0.88 


41.30 


-3.08 


-3.71 


early 






320 


0.14 


1 


0.05 


gal 


23.11 


1.04 


39.40 


-3.53 


-3.25 


late 


E (0.14) 


0.14 b 


323 


0.51 


1 


6.48 


agn2 


11.69 


3.65 


42.82 


-1.37 


-1.55 


late 




326 


0.36 


2 


0.05 








40.29 


-3.00 


-3.24 


late 






332 


0.56 


1 


0.04 


gal 


16.42 


0.71 


40.70 


-3.13 


-3.54 


late 






333 


0.38 


1 


0.86 


agn2 


0.95 


0.97 


41.62 


-1.96 


-2.01 


late 




0.38 b 


337 


0.90 


1 


0.03 


gal 


12.22 


0.90 


41.04 


-3.31 


-3.79 


late 






339 


0.25 


1 


0.05 


gal 


22.38 


1.14 


39.95 


-3.82 


-3.84 


late 


E (0.25) 


0.25 b 


344 


0.09 


4 


15.70 


agnl 


4.26 


0.99 


41.50 


-1.06 


-0.78 


ir 




346 


1.02 


2 


0.02 


gal 


7.37 


0.69 


40.98 


-2.80 


-3.76 


late 




1.02 b 


349 


1.12 


4 


0.53 


agnl 


3.71 


1.23 


42.56 


-1.14 


-2.45 


early 




1.25 bt 


351 


0.94 


1 


0.09 


gal 


7.85 


0.50 


41.59 


-2.64 


-3.53 


late 




0.94 b 


352 


0.91 


4 


0.03 


gal 


17.62 


0.69 


41.13 


-3.12 


-3.84 


late 






353 


0.42 


1 


0.10 


gal 


11.41 


1.05 


40.80 


-2.97 


-3.06 


late 


C (0.42) 


0.42 b 


354 


0.57 


2 


0.09 








41.05 


-2.53 


-2.70 


early 




367 


0.78 


4 


0.17 


agn2 


2.27 


1.19 


41.68 


-1.65 


-2.41 


early 






373 


0.48 


1 


0.15 


agn2 


22.40 


0.52 


41.11 


-2.35 


-2.68 


early 






378 


1.08 


2 


0.04 








41.40 


-2.35 


-2.98 


late 




1.08 b 


383 


0.17 


4 


0.05 


gal 


25.18 


1.47 


39.63 


-4.05 


-3.88 


early 


A (0.11) 


0.11 b 



TABLE 11 — Continued 



XID 


z 


Flag 


Fx 


Class 


log(oddsl) 


log(odds2) 


log(Lx) 


log(F x /F R ) 


log(F x /F K ) 


Type 


OBXF ID 


ZG07 


384 


1.02 


2 


0.09 








41.70 


-2.30 


-3.39 


6arly 






387 


0.97 


4 


0.08 








41.59 


-2 20 


-2.99 






1.01 b t 


389 


0.56 


I 


0.05 








40.79 


-2.52 


-3.15 


late 








n An 

U.4U 


4 


n n^ 


gal 


14 47 




40 48 


-3 24 


-o.ou 




Hi ^U.4_L J 


41^ 


401 


0.94 


1 


0.07 


y'a 1 


7.52 


0.63 


41.52 


-2.63 


-3.44 


early 




404 


0.11 


1 


0.07 


gal 


25.62 


1.11 


39.27 


-3.76 


-3.56 


late 


A (0.11) 




405 


0.94 


4 


0.05 








41.38 


-2.80 


-3.56 


late 




414 


0.80 


1 


0.09 








41.45 


-2.46 


-2.91 


early 






418 


0.28 


1 


0.09 


gal 


17.04 


1.26 


40.34 


-3.50 


-3.47 


late 




0.28 b 


426 


0.16 


2 


0.08 


agn2 


12.68 


0.75 


39.75 


-2.17 


-1.80 


ir 






428 


0.30 


2 


0.10 


gal 


14.84 


1.05 


40.47 


-3.23 


-3.13 


late 


E (0.30) 




433 


1.04 


4 


0.08 


gal 


7.72 


0.69 


41.66 


-2.70 


-3.61 


early 




1.02 b 


437 


0.82 


4 


0.85 


agnl 


2.18 


0.54 


42.44 


-1.80 


-2.37 


late 




0.84 b 


A £3 
4DO 


U. ( o 


,1 
'i 


U.oo 


; ign2 


Z.l 1 


u.oy 


A 1 

41. yo 


9 1 1 
-Z. LO 


9 7/1 
-Z. /4 


early 




n QA b 


454 


0.46 


4 


0.39 








41.48 


-2.26 


-2.45 


early 




0.46 b 




u.u / 


I 

i 


U.Ol 








^7 


-O.UO 


9 7S 


late 


(n 07^ 




462 


0.51 


1 


0.03 








40.44 


-3.21 


-3.36 


early 






466 


0.44 


1 


0.09 


gal 


13.55 


1.16 


40.81 


-3.15 


-3.65 


early 






471 


1.17 


2 


0.10 


gal 


7.89 


0.64 


41.90 


-2.71 


-3.71 


early 






473 


0.31 


1 


24.90 


agnl 


7.71 


3.56 


42.88 


-1.12 


-1.25 


late 






477 


0.44 


1 


0.10 








40.86 


-2.67 


-2.82 


early 






478 


0.08 


i 


2.27 


agn2 


4.40 


2.84 


40.55 


-1.03 


-0.82 


early 






480 


0.46 


1 


0.20 








41.19 


-2.28 


-2.77 


early 







Note. — XID is ID number from Alexander et al. (2003). z gives the adopted redshift. Flag describes the source and quality of the redshift: (1 = high-quality spectroscopic redshift, 
2=unknown-quality spectroscopic redshift, 3=low-quality spectroscopic redshift, 4=high-quality photometric redshift, 5=low-quality photometric redshift). Fx gives the 0.5-2.0 keV 
X-ray flux from Alexander et al. in units of 10 — 15 erg cm -2 s — 1 . Class gives the results of the Bayesian classification, if any. log(oddsl) and log(odds2) give the logarithms of the 
Baycsian odds ratios for the alternate models (agnl and agn2 in the case of galaxies, galaxies and agn2 in the case of agnl, and agnl and galaxies in the case of agn2). Lx is the 0.5-2.0 
keV luminosity of the source in erg s _1 . log(Fx/Fji) and log(Fx / Fjf ) give t he X-ray/R-band and X-ray /K-band flux ratios. Type gives the SED type computed in the photometric 
analysis. OBXF ID gives the spectral type and redshift of the sources found in Hornschcmcicr et al. (2003). zgo7 gives the redshift from Georgakakis et al. (2007), if present. 
a source is in the G07 infrared-faint sample 
k source is in the G07 infrared-bright sample 
T the G07 redshift is a photometric redshift 



TABLE 12 
CDF-S X-ray Sample 



05 



XID 


z 


Flag 


Fx 


Class 


log(oddsl) 


log(odds2) 


log(Lx) 


log(F x /F R ) 


log(F x /F K ) 


Type 


29 


0.57 


1 


0.24 


gal 


7.12 


0.81 


41.51 


-3.25 




early 


44 


0.57 


1 


0.60 








41.91 


-2.96 




early 


53 


0.67 


1 


0.14 








41.43 


-2.57 


-2.95 


early 


60 


0.54 


1 


0.23 








41.44 


-2.53 




late 


73 


0.42 


1 


0.10 


gal 


13.41 


1.22 


40.78 


-3.62 


-3.81 


early 


75 


1.00 


1 


0.08 


gal 


6.74 


0.48 


41.61 


-2.81 


-3.45 


late 


80 


0.58 


1 


0.23 


agn2 


6.17 


0.77 


41.49 


-2.21 


-2.38 


early 


83 


0.68 


1 


0.06 


gal 


9.35 


0.76 


41.11 


-3.00 


-3.39 


early 


84 


1.03 


1 


0.15 








41.93 


-2.30 


-3.28 


late 


88 


0.60 


1 


2.48 


agn2 


21.15 


2.29 


42.58 


-1.87 


-2.12 


late 


94 


0.12 


1 


0.73 


agn2 


30.05 


2.26 


40.45 


-1.43 


-1.69 


late 


103 


0.68 


1 


1.76 


agnl 


3.64 


1.47 


42.54 


-1.68 


-2.06 


early 


106 


0.67 


1 


0.27 


agn2 


1.67 


0.50 


41.72 


-1.98 


-2.45 


early 


113 


0.52 


1 


0.09 








40.97 


-2.45 


-2.49 


early 


115 


0.34 


1 


0.06 


gal 


18.63 


1.06 


40.33 


-3.66 


-3.84 


early 


117 


0.57 


1 


4.50 


agnl 


7.26 


3.69 


42.77 


-1.19 


-1.46 


late 


118 


1.10 


1 


0.55 


agn2 


38.53 


1.68 


42.56 


-1.80 


-2.60 


early 


121 


0.73 


1 


0.28 


gal 


6.78 


0.80 


41.84 


-2.97 


-3.38 


early 


122 


0.18 


1 


0.36 


agn2 


5.20 


2.97 


40.52 


-0.91 


-0.98 


ir 


124 


0.96 


1 


0.08 








41.56 


-2.51 


-3.06 


late 


126 


0.59 


4 


0.33 


agn2 


5.84 


2.85 


41.67 


-1.17 


-1.55 


early 


129 


0.23 


1 


0.22 


gal 


14.37 


1.15 


40.54 


-3.23 


-3.19 


early 


134 


0.60 


1 


0.07 








41.00 


-2.64 


-3.13 


early 


146 


0.73 


1 


0.09 








41.33 


-2.44 


-3.26 


early 


149 


0.13 


1 


0.23 








40.02 


-2.73 


-2.50 


late 


152 


0.25 


1 


0.06 


gal 


15.29 


0.90 


40.06 


-3.36 


-3.08 


ir 


155 


0.53 


1 


1.85 


agnl 


3.26 


1.37 


42.32 


-1.82 


-1.88 


early 


158 


0.36 


1 


0.09 


gal 


11.29 


0.73 


40.61 


-2.94 


-2.92 


late 


159 


1.05 


1 


0.91 


agn2 


4.15 


3.43 


42.73 


-1.00 


-2.01 


late 


161 


0.73 


1 


0.31 








41.87 


-3.09 




late 


162 


1.04 


1 


0.08 


agn2 


8.43 


0.78 


41.65 


-2.17 


-2.54 


late 


167 


0.58 


1 


0.09 


gal 


11.41 


1.00 


41.06 


-3.28 


-3.51 


early 


169 


1.02 


4 


0.05 








41.43 


-3.65 




late 


171 


1.03 


1 


0.18 


agn2 


0.64 


1.56 


42.00 


-1.83 


-1.72 


late 


176 


0.73 


1 


0.61 


agn2 


6.11 


1.71 


42.17 


-1.79 


-2.31 


early 


179 


0.66 


1 


0.49 


agn2 


30.04 


1.77 


41.97 


-1.87 


-2.13 


late 


181 


0.74 


1 


0.25 








41.79 


-2.32 


-2.61 


early 


182 


0.21 


1 


0.67 


gal 


10.83 


1.08 


40.95 


-3.13 


-3.06 


early 


184 


0.62 


1 


1.32 








42.33 


-2.86 




early 


185 


0.18 


4 


0.06 


agn2 


11.38 


0.60 


39.74 


-2.33 


-1.98 


ir 


189 


0.08 


1 


0.81 


gal 


21.22 


2.08 


40.06 


-3.60 


-3.42 


late 


190 


1.02 


1 


0.06 








41.48 


-2.09 


-2.70 


late 


192 


0.08 


1 


0.57 


gal 


21.21 


1.53 


39.91 


-3.56 


-3.43 


late 


193 


0.96 


1 


1.81 


agnl 


6.23 


3.07 


42.93 


-1.27 


-1.96 


late 


196 


0.67 


1 


0.07 








41.13 


-2.74 


-2.75 


late 


203 


1.18 


4 


0.04 


gal 


6.08 


0.72 


41.54 


-3.29 


-4.12 


early 


907 


o i n 

U. 1U 


i 
i 


o i s 


gal 


97 ^o 


1 9^ 

i.zo 


oy.oy 


A OQ 


A 09 


early 


210 


0.83 


1 


0.07 








41.38 


-2.18 


-2.71 


early 


212 


0.67 


1 


0.17 








41.53 


-2.52 


-2.80 


early 


214 


0.84 


1 


1.51 


agnl 


5.35 


1.98 


42.71 


-1.22 


-1.63 


late 


224 


0.55 


1 


0.12 


gal 


9.46 


0.88 


41.14 


-2.90 


-3.31 


late 


225 


0.55 


1 


0.06 


gal 


11.56 


1.05 


40.83 


-3.21 


-3.62 


late 


226 


1.00 


1 


0.12 








41.79 


-2.19 


-2.71 


late 


227 


0.67 


1 


4.01 


agnl 


7.49 


2.95 


42.89 


-0.90 


-1.19 


early 


229 


0.67 


1 


1.86 


agnl 


3.40 


1.91 


42.55 


-1.59 




late 



TABLE 12 — Continued 



XID 


z 


Flag 


Fx 


Class 


log(oddsl) 


log(odds2) 


log(Lx) 


log(F x /F R ) 


log(F x /F K ) 


Type 


236 


0.46 


1 


u. 


.UK 


gal 


9.70 


0.59 


40.78 


-2.84 


-2.91 


ir 


238 


0.24 


1 


0. 


.06 


gal 


17.46 


1.25 


40.05 


-3.42 


-3.66 


late 


244 


0.58 


1 


u. 


1 9 








41.21 


-2.54 


-2.99 


early 


247 


0.62 


2 


U. 


.0 1 


agn2 


22.40 


0.76 


41.95 


0.59 




early 


248 


0.67 


1 


U- 


,Uo 


gal 


10.45 


1.02 


40.83 


-3.35 


-3.74 


early 


256 


0.74 


1 


Q 
O. 


.on 


agnl 


9.56 


4.09 


42.96 


-0.55 


-0.88 


ir 


260 


0.27 


1 


U. 


.01 


gal 


7.85 


0.95 


41.08 


-2.68 


-2.82 


early 


262 


0.42 


1 


U- 


. lo 


gal 


10.01 


0.71 


40.89 


-2.91 


-3.05 


late 


263 


1.01 


1 


0. 


.56 


agn2 


12.02 


2.67 


42.48 


-1.53 


-1.85 


early 


265 


0.46 


1 


U. 




gal 


9.69 


0.88 


41.04 


-3.48 




late 


266 


0.44 


1 


0. 


.09 


gal 


10.56 


0.81 


40.80 


-2.93 


-3.13 


late 


267 


0.10 


1 


0. 


.11 








39.45 


-2.96 




ir 




U.DD 


i 

i 


0. 


.22 










-o.oz 




early 


271 


1.18 


1 


0. 


.12 








41.96 


-2.32 


-3.04 


late 


273 


1.11 


1 


0. 


.09 








41.79 


-2.30 


-3.02 


late 


276 


0.67 


1 


0. 


.15 








41.47 


-2.99 


-3.43 


early 


277 


1.15 


4 


0. 


.78 


agnl 


2.04 


1.16 


42.76 


-0.31 




late 


292 


0.37 


1 


0. 


.22 








41.01 


-2.70 


-2.42 


ir 


293 


1.14 


1 


0. 


.18 


gal 


8.46 


0.49 


42.13 


-4.43 


-2.47 


early 


303 


0.25 


1 


2. 


.70 








41.71 


-1.95 




early 



o 

o 

Note. — Columns are as in Table [Tl] CO 

o 



